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Abstract 

We initiate the study of noncharacteristic boundary layers in hyperbolic-parabolic 
problems with Neumann boundary conditions. More generally, we study boundary lay- 
ers with mixed Dirichlet-Neumann boundary conditions where the number of Dirichlet 
conditions is fewer than the number of hyperbolic characteristic modes entering the do- 
main, that is, the number of boundary conditions needed to specify an outer hyperbolic 
solution. We have shown previously that this situation prevents the usual WKB ap- 
proximation involving an outer solution with pure Dirichlet conditions. It also rules out 
the usual maximal estimates for the linearization of the hyperbolic-parabolic problem 
about the boundary layer. 

Here we show that for linear, constant-coefficient, hyperbolic-parabolic problems one 
obtains a reduced hyperbolic problem satisfying Neumann or mixed Dirichlet-Neumann 
rather than Dirichlet boundary conditions. When this hyperbolic problem can be solved, 
a unique formal boundary-layer expansion can be constructed. In the extreme case 
of pure Neumann conditions and totally incoming characteristics, we carry out a full 
analysis of the quasilinear case, obtaining a boundary-layer approximation to all orders 
with a rigorous error analysis. As a corollary we characterize the small viscosity limit for 
this problem. The analysis shows that although the associated linearized hyperbolic and 
hyperbolic-parabolic problems do not satisfy the usual maximal estimates for Dirichlet 
conditions, they do satisfy analogous versions with losses. 
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1 Introduction 

In the study of noncharacteristic boundary layers of hyperbolic-parabolic systems, phys- 
ical applications motivate the inclusion of Neumann boundary conditions along with the 
usual Dirichlet boundary conditions that have traditionally been considered for such prob- 
lems (see, e.g., [GSl [R2, R3] and rererences therein). In particular, as discussed in [NZ1, 
INZ2j IGMWZ5| [R] , suction- induced drag reduction along an airfoifj] is typically modeled by 
the compressible Navier-Stokes equations 

' dtp + div(pu) = 
dt(pu) + div(pu*u) + Vp = epAu + e(p + r/)Vdivn 
d t (pE) + dw((pE + p)u) = kAT + e/idiv((u • V)u) 

+ e(p + rj)V(u ■ divu) 

1 See [H[Br], or NASA site |http://www.dfrc.nasa.gov/Gallery/photo/F-16XL2/index.html| 
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on an exterior domain f2, with no-slip suction-type boundary conditions on the velocity, 
urldn = 0, u u \qq = V(x) < 0, and either prescribed or insulative boundary conditions on 
the temperature, T\ m = T waU {x) or d u T\ dn = 0. 

The study of such mixed-type boundary layer problems was initiated in [GMWZ5, 
GMWZ6 for certain combinations of Dirichlet and Neumann boundary conditions in the 
viscous problem. However, the ansatz used there, which assumes that the residual hyper- 
bolic problem should have only Dirichlet boundary conditions, breaks down when there 
are too many Neumann conditions in the viscous problem - more precisely, when there 
are too few Dirichlet conditions, in the sense that the number of scalar Dirichlet condi- 
tions in the viscous problem is strictly less than the "correct" number of residual boundary 
conditions for the hyperbolic problem. In such cases, the construction in [GMWZ5] of 
"C-manifolds" of reachable states determining Dirichlet boundary conditions for the outer, 
hyperbolic solution fails, due to a lack of transversality, as a consequence of which (together 
with the low- frequency decomposition of [R2j ) the maximal linearized estimates used in 
[GMWZ5, GMWZ6] to establish rigorous convergence may be shown to fail as well. As 
noted in [R], the case of (II. ip with incoming supersonic velocity falls into this category, so 
is not accessible by the techniques developed up to now. 

Clearly, in such cases, a new analysis is required. Several questions arise, including: 

(1) Does the hyperbolic-parabolic problem have a solution on a fixed time interval 
independent of e? 

(2) Is there a residual hyperbolic problem whose solution gives the small viscosity limit of 
solutions to the hyperbolic-parabolic problem? In particular, what are the correct residual 
hyperbolic boundary conditions? And, are these uniquely determined? 

(3) What are the maximal linearized estimates that we may expect in this context, both 
for the residual hyperbolic and full hyperbolic-parabolic problem? 

In this paper, we answer these questions completely in the extreme case of pure Neumann 
boundary conditions and totally incoming hyperbolic characteristic modes, showing that 
there is a reduced hyperbolic problem with Neumann instead of Dirichlet conditions, and 
that in place of the standard Dirichlet-type linearized estimates for the reduced hyperbolic 
and full hyperbolic-parabolic systems, there hold modified versions with losses, sufficient to 
close a rigorous convergence argument. As a corollary we characterize the small viscosity 
limit for the quasilinear problem. 

In the general, linear constant-coefficient case, we present two approaches to constructing 
a formal boundary-layer expansion to all orders of the solution to the hyperbolic-parabolic 
problem. In general the reduced hyperbolic (outer) problem features mixed Dirichlet- 
Neumann boundary conditions. In the pure Neumann case we prove that the exact and 
approximate solutions to the hyperbolic-parabolic problem are close when e is small. 

Our results motivate the further study of first-order hyperbolic initial-boundary-value 
problems with Neumann or mixed Neumann-Dirichlet boundary conditions. This is at first 
sight a counterintuitive problem, since the normal derivative on the boundary is not con- 
trolled by the usual hyperbolic solution theory, and it does not seem to have received much 
attention before now. We regard this as one of the most interesting aspects of the analysis. 
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1.1 Linear systems with Neumann boundary conditions 

First we examine a linear problem for which the above questions have a positive, and 
rather simple, answer. Let us consider the parabolic boundary value problem on M.^ 1 := 



{x = (x',x d ) = (x ,x",x d ) £ R d+1 : x d > 0}: 

Lu = f + eA x u in {x d > 0}, (1.2) 
9dU\ Xd =o = 0, (1.3) 
u\t<o = 0, (1.4) 

where L is a symmetric hyperbolic operator with constant coefficients 

d 

L = d t + y^Ajdj, t = x 

3=1 

and / £ H°°(M^~ d ) with f\t<o = 0- The N x N matrices Aj are constant (for now), and 
the boundary is noncharacteristic: 

det A d / 0. 

We look for an approximate solution of the form 

u e (x) = u (x) + eu\(x, -^) + e 2 u 2 (x, -^) + . . . 



Uj(x,z) = uAx) + u*(x',z), j>l, 



with the usual profiles 

where Uj is an "outer" solution, and u*j is a boundary layer profile which goes to as z — > oo. 

Remark 1.1. One could postulate a more general profile uo(x,z) = u (x) + Uq(x',z) at 
level j = 0; however, the resulting e _1 order profile equations A d d z u^ — B^Uq = 0, with 
boundary condition 9 z («q)| z=0 = would give then 8 z Uq = 0, recovering the assumption 
u = u (x). 

The profile equation obtained at the order e° is 

Luq + A d d z ui -d 2 z ui = f. 
which leads to the two equations for uq and u\: 

Lu = f (1.5) 

and 

A d d z u\ - d 2 z u\ = 0. (1.6) 
The boundary condition (j 1 . 3 [) gives at the order e°: 

(9dUo)\x d =o + {d z u\)\ z= Q = 0. 
Hence the solution to the boundary layer equation fjl .6[) is 

uUx, z) = -e zAd A- l d d u (x\ 0). (1.7) 
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It follows that u[ is decreasing at +00 if and only if d d uo\ Xd= o lies in E_(^4d)> the negative 
eigenspace of A d : 

d d u 0{xd=Q eE_(A d ). (1.8) 

But uo satisfies Luq = /; thus 

d-l 

d d u = -A2 l Y, A Au + A d 1 f 


and the condition (jl.8p is equivalent to 

Hu 0]Xd=0 G AjVU d =o + M40, (1-9) 

where If is the tangential operator iJ := A^" 1 Ylo ^ Ajdj. So we are led to solve the mixed 
problem 

Lu = fm{x d >0}, (1.10) 

Hu olxd=0 G AjV|*a=o+MAi). (1-11) 
n | 4<0 = 0. (1.12) 

(The boundary conditions may be rephrased via projections as described in Remark 11.41 ) 

To solve this problem introduce the unknown v := Huo, which is the solution of the 
symmetric hyperbolic problem with dissipative boundary conditions 

Hv + 8 d v = H{A- 1 f)m{x d >0}, (1.13) 
v\ Xd =o G A^f lxd=0 +E^(A d ), (1.14) 
v [t<0 = 0. (1.15) 

Hence v is completely determined; thus uq is also uniquely determined as the unique solution 
of 

Hu = v, -u |i<o = 

(here considered as an initial- value problem defined on slices x d = constant). Then u\ is 
uniquely determined by formula (jl.7p . and decays to zero at +00. 

The construction follows the same pattern for the next terms. For example, setting 
L' = dt + Ajdj we obtain at the order e 1 the profile equation 

Lu x + L'u\ + A d d z u 2 - d 2 z u 2 = Au . 

which leads to the two equations for ^ and u\: 

Lu x = Au (1.16) 

and 

A d d z u* 2 - d 2 z u* 2 = -L'u\. (1.17) 
The boundary condition fjl .3|) gives at the order e : 

(9dUi)| Xd =o + {d z u* 2 )\ z=0 = 0. 
One can solve as before these equations which gives a unique solution for u x and u\. 
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Theorem 1.2. u £ (x) = uq{x) +eu\(x, Xd/e) H Ufc(aJ, s^/e) + 0(e fc ) m L 2 ({— oo,T] x 

M+) /or all given T > and a// fc G N as e -> 0. 

Proof. Since we can construct an approximate solution to any order, it is sufficient to prove 
an estimate of |n|^2(n T ), where ft = (— oo,T] x M.+ , for the solution u to the problem 
(jl.2p(jl.3p(|l,4p . First we estimate the normal derivative. Applying dd to the equation 
(|1.2p and using condition (11. 3p leads to a hyperbolic-parabolic problem with a homogenous 
Dirichlet boundary condition for ddU. A simple integration by parts yields (with |«L = 
|e" 7 *u| L 2 (nr) ): 

elVxdd^ + 7|9du|^ < 7~ 1 |5d/lr 

Going back to the system (jl.2p . taking the product on the left by u, and integrating by 
parts leads to 

e\Vu\ 2 + 7|n|^ < t" 1 !/!^ + |it| 7 |5 d u| 7 . 
Hence using the previous estimate one gets 

e |Vu| 7 + 7 M 7 <T l \f\ 2 , + T z \d d f\ 2 v 

and finally 

M 7 ^7 _1 |/l7 + 7" 2 |^/l 7 - (1-18) 

Applying the estimate (|1.18p to the error w = u £ — u £ approx , with the function / replaced 
by 0(e r ) and ddf replaced by 0{e r ~ 1 ) for r chosen large enough (i.e., r > 2), proves the 
theorem. □ 

An analogous result with convergence in L 2 replaced by convergence in L°° can easily 
be obtained after getting higher derivative estimates. 

Remark 1.3. The approach followed here is similar to the idea of "filtering" introduced 
by Serre [Selj in the somewhat different context of second-order hyperbolic problems with 
variational structurejl in which a degenerate problem is decomposed into the composition 
of problems of standard type, each inducing its own losses/gains. 

1.2 Quasilinear systems with Neumann boundary conditions 

Next we derive a candidate for the residual hyperbolic problem in the quasilinear case. 
Consider the nonlinear parabolic problem 

L u (u) = f + eA x u in {x d > 0}, (1.19) 
ddU\ Xd=0 = 0, (1.20) 
u\ t<0 = 0. (1.21) 

where L u is a symmetric hyperbolic operator 

d 



L u = d t + ^2Aj(u)dj, 



Also featuring Neumann, or "free," boundary conditions. 
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and / G H°° (M. 1 ^~ d ) with /|t<o = 0. The matrices Aj are smooth and symmetric, and the 
boundary is noncharacteristic: 

detA d (u)^0, G R N . 

Again we expect an expansion of the form 

u £ (x) = u (x) + eui(x, ^) + e 2 u 2 (x, ^) + . . . , 

that is, a "weak" layer of order e in amplitude. This may be deduced exactly as in the linear 
constant-coefficient case, by examination of the order e _1 profile equations as described in 
Remark 11.11 

The equations for the terms of order e° give 

L Uo u + A d {u G )d z u\ - d 2 z u\ = f. 

This equation splits into two parts 

A d (u olxd=Q )d z ul - dlu\ = (1.22) 

and 

L Uo u = f, (1.23) 
and the boundary condition at the order s° is still 

(d d u ) lxd=0 + (d z ul) lz=0 . (1.24) 

The solution to the boundary layer equation f)1.22|) is 

u\{x',z) = -e^^'^A-^uoix'^^ddUoix^O). (1.25) 

This solution goes to at +oo if and only if 

d d u (x',0) G E_(Ai(u (x',0))). (1-26) 

Using the equation fjl .231) we rewrite this condition: 

H uo (u ) G A d 1 (u (x , ,0))f lxd=0 + E_(A d (u (x , ,0))). 

with H u := Adiu)' 1 L u — d d - Writing instead 

d-l 

L' u :=d t + Y J A ] {u)d r 
l 

we obtain the following hyperbolic boundary problem obtained for uq: 

L u {u) = f in(-oo,T] xR d + (1.27) 
L'uiu) G /|, d=0 +E_(A d (n)) on{x d = 0}, (1.28) 
u lt<0 = 0. (1.29) 
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Remark 1.4. We do not know if this problem is well-posed in general. The boundary 
conditions (|1.28p are unusual; they can be as rephrased as 

ir + (A d (u)) (L' u (u) - f\ Xd=Q ) = on {x d = 0} or 
ir + (A d (u))d d u = on {x d = 0} 

(equivalently, 7r + d d u = 0), where tt + (A(u)) is the projection onto E + (A d {u)) along E^(A d (u)). 
Yet, in the constant coefficient linear case the corresponding problem (ll.lOp . (II. lip . (|1.12p 
turns out to have a unique natural solution. 

In the totally incoming case where A d (u) > and thus E^(A d (u)) = 0, one can solve 
(jl.30p by first solving a hyperbolic system on the boundary, as we describe further be- 
low. A high-order approximate solution to the hyperbolic-parabolic problem (|1.19p can be 
constructed, and the small viscosity limit can be completely analyzed. 

1.3 Assumptions and main result. 

Our main result treats a quasilinear hyperbolic-parabolic problem where the questions 
posed at the beginning can be answered completely, the case where all characteristics for 
the hyperbolic problem are incoming: A d (u) > 0. We study the forward problem on 

if 1 : = i x = ( x '> x d) = ( x 0, x ", x d) e R d+1 : x d > 0}: 

d 

£{u e ) := ^2,Aj(u)d Xj u - eAu = f 

i=° (1.31) 
d Xd u\x d =o = 
u = in xo < 

where the Aj are N x N matrices (not necessarily symmetric), A d (u) > 0, and Aq = I. 
The approximate solution, which is constructed in section (27TJ has the form 

u a £ {x) =u°(x)+eu 1 {x) + --- + e M u M {x) (1.32) 

and satisfies 

d 

8{u a ) := A j (u a )d x .u a - eAu a = f + e M R E 

j=0 (1.33) 
d Xd u a \x d =o = 
u a = in x < 0. 

As a consequence of the totally incoming assumption, there is no fast transition layer in 
u a . Nevertheless, the nonlinear stability of u a and the analysis of the small viscosity limit 
turn out to be delicate questions, because the Evans function for this problem vanishes at 
zero frequency. Thus, u a can be expected to be at best "weakly stable" . 

The low frequency Evans function is computed explicitly in section [2~4l and its degeneracy 
near is precisely estimated H This estimate allows us to construct degenerate Kreiss 

3 Outside a neighborhood of zero frequency, the Evans function is nonvanishing by (JT7T5J ; recall that the 
layer in the totally incoming case is constant, so the analysis of Section fl - II applies. 
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symmetrizers at the symbol level in section 12.51 and these symmetrizers are used there to 
prove resolvent estimates for the frozen coefficient linearized problemQ 

The resolvent estimates are quantized in section 12.61 using the pseudodifferential calculi 
outlined in the Appendix. This section provides the main variable coefficient L 2 estimate, 
Theorem 12.61 for the problem obtained by linearizing the original system (|1.31|) around the 
approximate solution u a . Fortunately, the I? estimate exhibits no loss of derivatives, but 
there is a loss of a factor of yfe when the boundary datum g = 0. This loss in the main 
estimate, which reflects the degeneracy in the Evans function, is the source of most of the 
technical difficulties in the paper, because it prevents us from absorbing terms that would 
otherwise be absorbed easily as "error terms" in the estimates. 

Higher derivative estimates are proved in section 12.71 using an appropriate enlarged 
system, and these estimates are then used in section [2751 to solve the nonlinear error equation 
satisfied by u £ — u a by Picard iteration. 

We let Qt '■= {x = (x',Xd) = (xQ,x",Xd) G M. d+l : Xd > 0, xq < T} and sometimes 
write t = xq. 

Assumption 1.5. I.) The N x N matrices Aj{u) in the system (|l,3ip are C°° and sym- 
metric, Aq = I, and Ad{u) > 0. Thus, in particular the boundary is noncharacteristic. 

II.) Let f G fP(IR + +1 ) for s large (as in Theorem \2.18\) . f = in t < 0, and let 
uo(x) G Qt denote the solution to the residual hyperbolic problem: 

d 

dtu + Aj(u )dju = f in x d > 

j=1 (1.34) 

ddUo\ Xd =o = 

uq = in t < 0. 

Assume that for x G £It the function uq takes values in a neighborhood of 0, U, such 
that for u G W, the hyperbolic operator dt + Ylj=i Aj(u)dj has semisimple characteristics of 
constant multiplicity. 

Remark 1.6. The positivity of Ad implies that the boundary condition in (jl.34p agrees with 
(|1.30p . Assumption II is a familiar condition implying that the hyperbolic system satisfies 
the "block structure" condition first formulated by Kreiss [K] for constructing symmetrizers. 
We could replace Assumption II by other weaker assumptions that imply block structure. 
We could also require that such an assumption holds only for x near Xd = with only minor 
changes in the proofs. 

Theorem 1.7. Under Assumption \1.5\ there exists an such that for < e < €q the 
parabolic problem (|l,3ip has an exact solution u £ on $7t of the form 

u e {x) = u a £ +e L v £ , (1.35) 

where u £ has the expansion (11.32P in which the leading term is the solution uq to the residual 
hyperbolic problem ([1.34p . The exponent L can be chosen as large as desired provided the 

4 Degenerate symmetrizers were used also in GMWZ2 , but there the degeneracy occurred in the elliptic 
bloc (Sp in (I2.46|) '). rather than the hyperbolic block. 
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approximate solution is constructed with sufficiently many terms (M(L) ) and in that case 
we have: 

\d a (v E ,ed d v £ )\ L oo < 1 (1.36) 
for \a\ < L, < e < eo- Here d = (do, . . . , d c i-i)- 

This Theorem is an immediate corollary of the more precisely stated Theorem 12.181 
which is phrased in terms of U = (v,eddv). 

Corollary 1.8 (Small viscosity limits). Let u e be the solution to the hyperbolic-parabolic 
system (jl.31j) . u a e the approximate solution fjl.32H to that system, and uo the solution to the 
residual hyperbolic problem (|1.34j) . Then 

\u E -u a e \ L oo {(lTo) <Ce L 
\u £ - uo\L°°(n TQ ) < Ce. 

1.4 Mixed boundary conditions: toward a general theory 

We conclude with a discussion of the case of mixed Dirichlet-Neumann boundary con- 
ditions in the linear constant-coefficient case, making contact with the previous work of 
[GMWZ5J. Consider again a linear constant-coefficient boundary value problem 

Lu = f + eA x u in {x,} > 0}, 

for L as in section [TTTjj with mixed boundary conditions 

^2ddu\ Xd=0 = g 2 

satisfying 



(1.38) 



rankri + rankl^ = rank ( p 1 ) = N. (1.39) 



2 

Let us suppose now that f,gi, and g 2 vanish in t < and satisfy high-order corner com- 
patibility conditions at t = 0, Xd = 0. We seek u such that u = in t < 0. 
We seek a formal boundary-layer expansion 

u b (x) = u (x, —) + +eu\{x, —) + e u 2 {x, —) + ... 

with profiles 

Uj(x, z) = Uj(x) + u*(x', z), j > 0, 

where Uj is an "outer" solution, and u*j is a boundary layer profile which goes to as z — > oo. 

Denote by rankri =: T> the number of Dirichlet conditions, rankl^ =: J\f the number 
of Neumann conditions, dimE + =: I the number of incoming modes, and dimE_ =: O the 
number of outgoing modes, so that 

V+M = Z + = N. 



Evidently, we can extend as in Sec. ll.3l to the nonsymmetric case, at the expense of further assumptions. 
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Henceforth, we may (and do) take F\ to be a T> x N matrix and T2 to be an J\f x N matrix. 
We divide the analysis into two cases: 

(i) V > X, or, equivalently, M < 0, and 

(ii) V < X, or, equivalently, M > O. 

The first case is the one considered in |GMWZ5] , and treated for problem in [Rj . The 
second includes the case of Neumann boundary conditions treated here, and also the case 
of problem (II. ip left untreated in |R r As we shall see, they have quite different behavior. 
We will see that in case (i) the reduced boundary condition on u is derived as a solvability 
condition for obtaining Uq, while in case (ii) Uq = and the reduced boundary condition on 
Uq is derived as a solvability condition for obtaining u\. We begin by recalling, with some 
simplifications possible for this linear problem, the treatment of case (i) in [GMWZ5]. 

Case (i). The general solution of A<id z UQ — d^u^, which decays to as z — > 00, has the 
form 

u* (x',z) = e AdZ d(x') (1.40) 

where d £ E_(^) is arbitrary (here and henceforth we suppress x'). The e _1 order bound- 
ary condition r2<9 z Uo(0) = implies 

d z u* Q (0) e ker (r 2 | E _ (Ad) ) and thus 1*5(0) G A" 1 ker {T 2 \ E _ (Ad) ) . (1.41) 

We make the following transversality assumption: 

(a) To has full rank, namely J\f, on E_(Ay) 

\ , x (1-42) 

(b) Fi has full rank on X := A d ker (r2|E_(A d )) • 

Since dimE_(Ad) = O, Assumption II ,42l f a) implies dim(ker (F2\^,_(A d ))) = O — M and thus 
(fL42]) (b) implies 

dimTiX = O-M. (1.43) 

Since the subspace F\X C MP and T> = X + O — N, F±X is equal to the null space of some 
X x V matrix, call it /C. Now use the order e° Dirichlet condition 

F 1 (u o (0) + u* Q (0))=9l (1-44) 
to see that there exists Uq(0) £ X satisfying (jl.44p if and only if 

ri(«o(o))-0i€rix (1.45) 

In other words 

f i(«o(0)) = gi, (1.46) 

where F% = JCF\ and gi = K,g\. Observe that F\ is an X x N matrix of rank X as required. 
The reduced hyperbolic problem is therefore 

Lu = f 

?iu.o = gi on x d = (1.47) 
u = in t < 0, 
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which is well-posed provided that the usual Kreiss Lopatinski conditiorjl is satisfied. Con- 
tinuing this process, one obtains an expansion to all orders. In this case, boundary layers 
are amplitude 0(1) and the reduced boundary conditions are purely Dirichlet. 
Case (ii). We now turn to case (ii), where we make the assumption 

T 2 is full rank on E_ (A d ). (1.48) 

Since M = rankl^ > O = dimE_(A^)), we find from the e _1 order profile equation 
T 2 d d u* (0) = 0, and the fact by (fOO]) that d d u* G E_(A d ), that 

d d u* = u* = 0. (1.49) 

Thus, the boundary-layer expansion features a weak layer of amplitude 0(e), just as in the 
full Neumann boundary condition case. This implies by the order e° boundary condition 
Fi^o = <7i, and the weak layer property uq = u , that the Dirichlet condition is inherited 
unchanged by the outer solution, as 

^im\x d =o = 91- (1-50) 
The order e° Neumann condition is 

r 2 (9dWoU d =o + d z ul\ z=0 ) = g 2 . (1-51) 

We deduce the reduced Neumann condition on uq as a solvability condition that allows us 
to find a solution d z u\\ z =o G E_(^4rf) of (|1.51|) . Recalling that T 2 is an M x N matrix, we 
denote by S the subspace of M"^ given by 

S = T 2 (E^(A d )). (1.52) 

By (|1.48p the dimension of S C R-^ is O. Thus, S coincides with the kernel of an (N—O) xM 
matrix. Choose one such matrix and call it M. 

By the definition of M, in order to find d z Ui\ z= o G K^(A d ) satisfying (|1.5ip we must 
have 

M (T 2 (d d u \ Xd = ) - g 2 ) = 0, (1.53) 

or in other words 

r 2 d d u \ Xd= o = g 2 , (1.54) 

where T 2 = MT 2 and g 2 = Mg 2 . As expected, T 2 is an (M — O) x N matrix of rank 
(N — O), giving us the remaining N — O boundary conditions needed (in addition to the 
D Dirichlet conditions) for the hyperbolic problem. 

Combining, we obtain the reduced hyperbolic boundary-value problem 

Lu = f in x d > 

(1.55) 

r2<9d«oU d =o = 92 

uq = in t < 0. 



6 In [GMWZ5] it is shown that both the Kreiss-Lopatinski and transversality conditions follow from a 
condition on the low-frequency behavior of an Evans function. 
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Remark 1.9. a) In the case of full Neumann boundary conditions we have M = N = I+O, 
and T2 is a nonsingular N x N matrix, which we may therefore always take to be In. Then 
we have S = E_ (Ad) (|1.52[) and we may take M = T 2 to be an (N — O) x N matrix whose 
rows span K + (Ad). 

b) In the totally incoming case with full Neumann boundary conditions we have O = 0, 
S = {0} C R N , and we can take M = Id N . So f 2 = T 2 = I. 

c) In the totally incoming case with one Neumann boundary condition, we have N = 1, 
V = N - 1, S = {0} C R 1 , and we may take M = 1. Thus, f 2 = T 2 , a 1 x N matrix. 

d) In the totally incoming case we have E_(Ad) = {0}; thus, our construction of the 
the approximate solution shows that u*-(x',z) = for all j. In other words, the layer is 
absent (or constant). 

e) In the situation T> = X on the boundary of case (i) , assuming (jl.48p , we find by the 
argument of case (ii) that the amplitude of boundary layers is 0(e). In other words, the 
layer is absent to lowest order also in this boundary case. 

By introducing variations on the method of Section 11.11 we discuss next two approaches 
to obtaining a well-posedness theory for problems of the form (11.55|) . When one has such a 
theory, one can proceed as in section 11.11 to construct the boundary layer expansion to any 
order. 

1.5 The reduced hyperbolic problem: approach based on Kreiss sym- 
metrizers. 

Substituting for ddUo the expression 

d-l 

d d u = -A- 1 (d t u Q +^ / A j d j u ) + A- 1 /, (1.56) 
3=1 

and taking the Laplace- Fourier transform with Laplace frequency 7 + ir, 7, r 6 R 1 , and 
Fourier frequency r] S R rf_1 , we convert the boundary operator appearing in (]1.56|) to the 
homogeneous degree one boundary symbol 

d-i 

-AdHl + ir + ^iVjAj) (1.57) 
3=1 

The matrix (7 + ir + Ylj=i ^VjAj), by symmetry of Aj, is invertible for 7 > with 0(7 _1 ) 
inverse. As we saw above T 2 is of full rank r := M — O; hence T' 2 ■= — f , 2^4^" 1 (7 + 
ir + ^2jZiiT]jAj) has the same rank for 7 > 0. Multiplying on the left by m(7, t, 77) : = 
(ir + 7 + \r]\)~ , we obtain a symbol homogeneous of degree zero 

f 2 (7,r,T?) := -m(7,T,7 ? )f 2 A- 1 ( 7 + iT + ^^ J 4 j )- (1-58) 

3^d 

The Neumann boundary conditions can be rewritten now as degree-zero Dirichlet conditions 
f 2(7, r, 77)m (7, t, t], 0) = G 2 (7, t, rj) 

h ~ ,~ \ (1.59) 

:= m(7,r,r/) (52(7, r, 77) - r 2 .4 d /(7,r,r/,0)J , 
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where " denotes Laplace-Fourier transform. 

With this rephrasing of the boundary conditions, the Laplace-Fourier transformed sys- 
tem becomes a hyperbolic boundary-value problem of the following form: 

d-i 

d d u + A d 1 (<y + it + ^2 iVj A j)uo = 1 /(7, t, rj, x d ) 

j=1 (1.60) 

rin (7,r,r/,0) = g x 

F2 (7, t, v)uo{l,T,7],0) = G2(7,t,t/) as in f|1.59|). 

Uniform estimates may be proved for (jl.60p using Kreiss symmetrizers (see, for example, 
\CP\ \BS\ IMet4| . and also Proposition IB. 1|) . provided that: (i) the boundary matrix T is 
uniformly well-conditioned, 

\T\,\T^\<C, (1.61) 
where is the pseudoinverse of V, and (ii) there holds the uniform Lopatinski condition: 



det ker 



fa^r^))'^ i A ^ + iT + i ^ A ^ >C>0 (1.62) 



for some C independent of (r, n) G 7 > 0. Here one defines the determinant by taking 
an orthonormal basis for each of the spaces appearing there. The condition thus expresses 
"uniform transversality" of those spaces for all such (7, r, rf). 

For discussion below, we recall also the weak Lopatinski condition, which is defined as 
in (|1.62|) . except that C 7 > is allowed to depend on 7 > 0. 

Assuming that the uniform Lopatinski condition is satisfied, we can use the following 
proposition to solve the outer hyperbolic problem. In the next proposition for 7 > 1 we let 

I/U.7 := l r )7,^r/(T-«7,^^d) r _, (1-63) 

L A {T,r),x d ) 

and we let {g)s,y denote the corresponding norm on the boundary. The block structure 
assumption made in the next proposition is satisfied by many of the important physical 
examples (see [MZ2| ): we shall omit further discussion of it here0 

Proposition 1.10. Suppose that L is an operator that can be conjugated to block structure 
in the sense of JMZ2J. Assuming well- conditioning (jl.6ip and uniform stability (|1.62p . there 
exist positive constants C , 70 and a unique solution of (jl.55p satisfying 

7l< 7 + « 7 < C (1/1^/7 + + (9i)l, + (92) 2 -i, 7 ) ■ (1-64) 

for 7 > 70. 



7 The block structure assumption can actually be avoided in the constant coefficient symmetric case by 
using the approach of [GMWZ8] . 
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Proof. For the problem (|1.60p with data (A d 1 f, gi, G2) one has the standard Kreiss estimate 
(ICBIBS]): 



l\u\l n + {u)l n < C ( J% + { 9l )l n + {G 2 )IA . (1.65) 



Existence for the problem (|1.6U|) follows from Proposition IB.ll which allows for pseudod- 
ifferential boundary conditions. The estimate (jl.64p now follows directly from (jl.65p and 
(|1.59p using \m\ ~ |r, 7,r/|~ 1 and 

(/U d =o)o, 7 <|/|o, 7 + |9, d /|o, 7 . (1.66) 

□ 

Assuming that the uniform Lopatinski condition is satisfied, we can solve the outer 
hyperbolic problem in this way and as in section 11.11 proceed to construct the boundary 
layer expansion to any order. The following proposition provides some information about 
when the weak and uniform Lopatinski conditions are satisfied by the problem (|1.60p . 

Lemma 1.11. Consider the problem (|1.60p . where the Aj are constant, real, symmetric 
N x N matrices. 

(a) In the totally incoming case with mixed boundary conditions or full Neumann bound- 
ary conditions, if the weak Lopatinski condition holds then the uniform Lopatinski condition 
holds. 

(b) Let d > 1. For full Neumann boundary conditions the weak Lopatinski condi- 
tion holds. The uniform Lopatinski condition can fail if the characteristics are not to- 
tally incoming. For example, it fails whenever there exists an eigenvalue uj{t — ij,rj) of 
— A^ir — 27 + Ylj=i AjTjj), analytic in r — ry, such that uj(t, rf) = and 8 t uj(t, rj) < for 
the chosen (t, rj). 

(c) For pure Neumann boundary conditions and d = 1 the uniform Lopatinski condition 
is satisfied. 

(d) In the totally incoming case with a single Neumann condition, the weak Lopatinski 
condition holds if and only if ( J (in this case a full N x N matrix) is invertible. 

V^ A d J 

(e) For mixed boundary conditions the weak Lopatinski condition holds only if ( ~ 

is full rank on E + (^4^). There are examples with mixed boundary conditions where weak 
Lopatinski fails and other examples where uniform Lopatinski holds. 

Proof, (a) In the totally incoming case E + ( J 4^" 1 (7 + it + ^X^j=i VjAj)) = C N . If the weak 
Lopatinski condition holds the determinant (jl.62p is ±1 for all 7 > 0. 

(b) In the full Neumann case T\ is absent and T2 is an (N — O) x N matrix whose 
rows span E_|_(yLj) (see Remark ll.9p . Since A^ l (^ + it + X)j=i ^Vj A j) is invertible for 7 > 
and E_|_( j 4^ 1 (7 + zt + ^^~! irjjAj)) an invariant subspace, we find that the weak Lopatinski 
condition is equivalent to T2 being full rank on E + (A ( ^ 1 (7 + it + rjjAj)) for 7 > 0. 
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Since the problem 



d-l 

d d w + AJ X ( 7 + it + im A i) w = „ 

j=i (1-67) 

f 2^U d =o = h 

is maximally dissipative, a simple energy estimate shows |ui(0)| < C\h\ when w G E + (A^ 1 (7+ 
ir + i Ylj=i VjAj)), so the full rank condition holds. 

In the case is not positive definite, the kernel space in fjl .62(1 must be nontrivial. 
Taking 7 = 0, \r)\ = 1 and choosing r from among the eigenvalues (77, 0) (here £<j = 0) of 
— X)i=i VjAj with corresponding eigenvector tj 7^ 0, we find that r' 2 (7, r, 77), or, equivalently, 
^2(7, T ) annihilates v. It can happen that v lies in the limit space as 7 — > of E + (^4^ 1 ( 7 + 
ir + i X]j=i VjAj))- The Cauchy-Riemann equations imply that this happens, for example, 
whenever there is an eigenvalue ui(t — 17, 77) of — ^ 1 (t— *7 + X)i=i Ajfjj), analytic in r — i 7 , 
such that w(r, 77) = and d T u(T, 77) < for the chosen (r, 77). 

Since u is also a limit of vectors lying in ker 1^2(7, r, 77) as 7 — > 0, we see that for such 
(t, 77) the determinant in (|1.62p converges to zero along some sequence j n — > 0. 

(c) When d = 1 and 7 > 0, we have ker 1^2(7, r) = kerf^A^" 1 and E + ( J 4^" 1 (7 + «r)) = 
E + (AA. Thus, both spaces are independent of (r, 7). The uniform Lopatinski condition 
now follows from the fact that I^ is full rank on E + (AA. 

(d) Recall from Remark ll.9l that T2 = T2 alxN matrix. The assertion follows by the ob- 
servation that in this case the real part of the determinant of f „ .1, . 1 . „ . s 

(e) The first assertion regarding mixed boundary conditions follows by inspection of the 
case 7 = 1, t = 0, 77 = 0. For the second assertion we refer to the examples given below. 

□ 

Remark 1.12. 1) When the problem (jl.60p only satisfies the weak Lopatinski condition, 
there is a still a chance of proving well-posedness for the reduced hyperbolic problem (|1 .55[) 
using degenerate Kreiss symmetrizers and constructing the WKB expansion. Indeed, several 
kinds of weakly stable problems have been studied successfully in this way (see, for example, 
. ICo2| ): typically the energy estimates exhibit a loss of derivatives. 

2) Glancing points are points (r, 77) where the matrix A^^ir + X^=i ^Vj^-j) has non- 
trivial Jordan blocks, or equivalently, where an eigenvalue Xj(£,r)) of YljZi r lj^ L j + £Ai> 
is stationary with respect to £. Such points always occur in d > 1, except in the totally 
incoming or totally outgoing cases, where they never occur (see [GMWZ6]). Example 11.11 
shows that the uniform Lopatinski condition can fail at glancing points. We know of no 
proof of well-posedness for the rescaled initial-boundary value problem in the case when the 
uniform Lopatinski condition fails in this way. (We present a different method in Appendix 
[D]for which this difficulty does not appear; see Example ID. 31 ) 

3) Example 11.31 shows that even weak stability can fail for the problem (11.60p . 



is 7 det 
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Example 1.1. Consider the simplest example of the first- order wave equation with drift a, 

A) l\ , (l + a 
A *={l oj' M= { -l + a 

with full Neumann boundary conditions, so that Y 2 = (l 0) . Then 
T' 2 = -t 2 A 2 \ 1 + ir + ir 1 A l ) = -(^ ^ 

which leads to the zero-order boundary matrix T 2 = — ir+ ^ + ^ ( i+a~ TT^J j ■ ^ Lemma 
\l.ll]( b) the weak Lopatinski condition holds. Applying the criterion of Lemma M.llV b). 
we find that the uniform Lopatinski condition fails at 77 = — 1, 7 = 0, r = 1, where 
lim 7 _> + ^+(^2 X (7 + ^ r + ^^1) = Span{(l, 1) T }. Moreover, the computation A±(£, 77) = 
a£ ± \/ S, 2 + rj 2 , where \± are the eigenvalues of ^A^+rjAi shows that d\±/d£ = a£ £ = 0, 
corresponding to failure at a glancing point, occurs only for a = /or t/iis choice of (r, 77) . 

Example 1.2. Next, consider the totally incoming problem 

Ai = K JV ,4 2 = Id, 
TOt/i mixed Dirichlet- Neumann conditions T% = ( * l), r 2 = r 2 = (l 0). Then, 

r A = ( Ti ] = ( Tu Ti2 

T'J \-t 2 A- 1 (-i + iT + ii 1 A 1 )) \-{l + iT) -irj 

is full rank on E+(A7 1 (7 + ir + mil)) = C 2 whenever 7^ det | ^ n ^\ 12 | = 7 + i(Y — 

?7rii), in particular for 7 > 0. Thus, we have weak Lopatinski stability of the zero-order 

boundary condition. ( - 1 J. By Lemma \l.ll\ the uniform Lopatinski condition also holds. 
\T 2 J 

Example 1.3. Finally, consider the totally incoming problem 

/0 1 a\ 
Ai = 1 1 , Ao = Id, 




uwi/i mixed Dirichlet- Neumann conditions F% = (l 1 6), r 2 = r 2 = ( n n I • Then, 




T'J ^-f 2 A 2 - 1 ( 7 + ir + Mi) 
is /it// ranA; on E + = C 2 w/ien its determinant is nonvanishing. An easy row reduction gives 
det(^)=det(( 7 + ir)Id-^^ « ) ) = 
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when (7 + ir)/irj is an eigenvalue of K ^\ , or (7 + ir)/irj = ab± ^ a2 J> 2 ± ^k _ Choosing 

a = 1, b = — 1, we obtain (7 + ir)Jir] = ~ 1 ^ v ^ ; or 7 + it = (—i =p \/3)(??/2), and i/ie 
Lopatinski condition is violated for 7 = ^f^/3r]/2, r = —rj/2. This shows that the weak 

Lopatinski condition can fail for the totally incoming case, even with ( ^}—\ J /wZZ ran/c. 



Example 1.4. 27ns Zasi example comes from a result by B. Fornet (see /FJj Z. ]F2j ). and 
shows that such types of Dirichlet- Neumann boundary conditions have a natural place in 
the theory of first order hyperbolic Cauchy problems with discontinuous coefficients. Let us 
consider the following scalar Cauchy problem in ID 

J d t u + a(x)d x u = f x £R,t>0 
\ u\ t =o = h 

where the real valued coefficient a(x) satisfies a(x) = —a < if x < and a{x) = (3 > if 
x > 0, with data h £ C°°(M), / G C°°(]R 2 ) compactly supported. This problem is of course 
not well-posed due to the lack of uniqueness. In order to select one solution, one can use 
for example a vanishing viscosity approach, and look for the limit of the solution u £ of 

J dtu + a(x)d x u — ed 2 x u = f x GR,t > 
\ u\ t =o = h 

as e — > 0. To study the convergence, the problem is written as an initial boundary value 
problem (or transmission problem) with u±(t,x) = u £ (t,±x) for x > and v £ = (u £ + ,u e _) T 
leading to the constant coefficient system 

d t v £ + Ad x v £ - ed 2 x v £ = (/+, /_) T in t > 0, x > (1.68) 

with 



A 



13 
a 



corresponding to totally incoming characteristic fields. The boundary conditions are Y\v = 
and T2d x v = on x = with 

r 1 = (i,-i), r 2 = (i,i). (i.69) 

The result is that v £ converges in L 2 ([0,T] x R + ) to the (unique) solution v° of the limit 
hyperbolic problem 

dt v + Ad x v° = (f + J^ T 
with the same boundary conditions 

r lV ^ x=0 = 0, r 2 (d x v\ =0 = 0, 

and initial conditions (h + ,h^) T . The fact that the problem is one dimensional helps a lot, 
and as a matter of fact, it is an example where the uniform Evans condition is satisfied (see 
\Fiy ). The convergence analysis also uses specific boundary layer expansions. One can find 
more general situations and examples in the paper lF0j with larger systems, still in ID. 
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Remark 1.13. Example 11.21 is an example of the mixed, totalling incoming case with one 
Neumann condition where the uniform Lopatinski condition holds. Recall that this case, 
corresponding to supersonic incoming flow with a Neumann condition on temperature, was 
left open in the study of boundary layers for the full compressible Euler equations (jl.ip in 
[R] . In Appendix ICl we provide a criterion (satisfied for example by ideal gases) for the 
uniform Lopatinski condition to be satisfied in that case. 

We point out that the well-conditioning of T, (jl.6ip . fails in many cases. In particular, 
in the totally incoming case, when there is even one Neumann condition, we find that T 
drops rank for 7 = at any values of r, ij for which r + Ylj^=d VjAj is not invertible, so that 
|rt| blows up as 7 — > 0. In this particular case, this may be remedied by simply multiplying 
r and data g both by to eliminate this difficulty at the expense of losses on the source; 
we explore this approach further in Remark 11.171 b elow . 

In this section, we have dealt entirely with construction of approximate solutions. Con- 
vergence to these solutions is a separate issue that requires estimates on the full hyperbolic- 
parabolic problem, estimates that we have for the moment only for the pure Neumann 
boundary, totally incoming case. This is an important direction for further investigation. 

1.6 Second approach based on solving a Cauchy problem on the bound- 
ary. 

We return now to the reduced hyperbolic problem in its original form (|l,55p in the 
general case of mixed-type boundary conditions, but assuming that we are in the totally 
incoming cas^l. Extensions to the general case are discussed in Appendix[Dl Writing u = uq, 
differentiating the Dirichlet boundary condition Tiu\ Xd= o = gi with respect to time, and 
making the usual substitution (|1.56p for d^uo, we obtain the boundary condition 

Bu\ Xrl= o = ( _ „ I , where B := ( ^\^-\ ] dt + ( ^ ?_i . ) d x . 

(1.70) 

The next proposition shows that sometimes this may be treated as a Cauchy problem in 
the tangential variables and solved for complete Dirichlet data u\ Xd= §. 

Definition 1.14. Let P (t, V ) := det (Yf^-i) r + £ f = \ (_f J-i^.) • We say 

that the system (]1.70p on the boundary is: 

a) evolutionary if the coefficient of dt is invertible. 

b) weakly hyperbolic if for any n £ the roots in r ofp(T,rj) = are real. 

Proposition 1.15. The system (jl.70p is both evolutionary and weakly hyperbolic if and 
only if the problem (|1.60p satisfies the weak Lopatinski condition. 

Proof. 1. First observe that the factor m in T 2 has no effect on the kernel space in (ll.62p . 
Suppose the weak Lopatinski condition holds. Taking r = 0, r] = 0, 7 > in (|l,62p . since 
the E+ space in (|1.62p is C^, we see that the coefficient of dt in (jl.70p is invertible. More 

8 Thus, this approach is relevant to the example of Rao discussed in Appendix (ICl) . 
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generally, the matrix ( - 1 ) is nonsingular when 7 > 0, and thus so is the matrix 
\r 2 (7,T,77)/ 



(-ivi- 1 ) (r " ll) + E ^ 





2. The argument can be reversed to prove the other direction. 



□ 



Weak hyperbolicity is not enough to guarantee well-posedness in H s spaces of the prob- 
lem (|1.70p . We refer to |BS] for a discussion of necessary and sufficient conditions for such 
well-posedness. One important sufficient condition for well-posedness is that the roots in r 
of p(r, 77) = are real and semisimple with constant multiplicities for r\ 7^ 0. This condition 
is verified for the system (jl.70p arising in the Rao example in Appendix [Cl 

In problems where the trace uq\ Xd =o = h can be found by solving (|1.7U|) . we can obtain 
the solution to the reduced hyperbolic problem (|I,55p by solving 

Lu = f in x d > 0, uo\ Xd= o = h, u = in t < 0. (1-71) 

This problem is maximally dissipative in the totally incoming case. 

We record the resulting bounds, which are to be compared to those of (|I.64p . 

Proposition 1.16. Suppose that L is an operator that can be conjugated to block structure 
in the sense of JMZ2 .J. Assuming that the roots in r of p{r,rj) = are real and semisimple 
with constant multiplicities for 77 7^ 0, there exist positive constants C , 70 and a unique 
solution of (|1.55p satisfying 

7H8 l7 + Ho, 7 < C (|/|o 2 7 /7 + \d.J\l.h 2 + (dt9i) 2 0r( h 2 + (h) 2 0r( h 2 ) ■ (1-72) 
for 7 > 70. 

Proof. Estimating the (-)o,7 norm of the trace of / at = 0, and using this to bound the 
f, A -\ , } in (|1.70p . we obtain from standard hyperbolic Cauchy estimates the 

92 ~L2A d fj 

bound (u| Xd=0 )§ (7 < C (|/|g i7 /7 2 + \d.J\l.ll 2 + {dt9i)l a h 2 ) + (hfo^h 2 ) , from which 
(|1.72p then follows by standard boundary value estimates for maximally dissipative systems. 

□ 

Remark 1.17. The bounds (jl.64[) obtained by method one in Proposition II . IQl are stronger 
than those of (ll.72p by factor 7/I7, r, rj\ in boundary terms gj and the term d Xd f coming 
from the trace of /. This reflects the well-conditioning hypothesis (jl.6ip made in Proposition 
11.101 but not in our derivation of (ll.72|) . Indeed, when well-conditioning fails (but the other 
hypotheses of Proposition 1 1 . 101 hold) one can apply method one to derive the bounds (|1.72p 
provided one can find for 7 > and r := rankTi + rankr2 an r x r matrix multiplier 
|m(7,r, 77)1 < C/7, such that the rescaled boundary condition 

( (7 + *r)ri \ 

m \-^A- d \ 1 + l r + Y:^A 3 )) 

satisfies the well-conditioning hypothesis (ll.6ip needed to obtain standard Kreiss-type 
bounds for the resulting rescaled boundary-value problem. One may check that this yields 
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exactly the bounds (jl.72p . Thus, this modification allows somewhat wider application 
of method one. For example, in the case of totally incoming characteristics, the uni- 
form Lopatinski condition is trivially satisfied, but (|1.61|) fails for the multiplier I7, r, 77I -1 , 
whereas for the multiplier m := the rescaled boundary condition mT = Id trivially 
satisfies (|l,6ip . and in favorable cases satisfies \m\ = | r 1 1 < C/7. Indeed, this can be 
recognized as the solution operator of the Cauchy problem on the boundary just described 
in method two. 

1.7 Discussion and open problems 

To summarize, following up on the analyses initiated in [GMWZ51 IGMWZ6] to accom- 
modate mixed Neumann-Dirichlet boundary conditions in the general theory of hyperbolic- 
parabolic boundary layers, we here investigate the case left open in those works that the 
number of incoming modes exceeds the number of Dirichlet conditions imposed on the full 
hyperbolic-parabolic solution. In this case, we find that (i) the resulting reduced, hyper- 
bolic, "outer problem" satisfies Neumann or mixed Neumann-Dirichlet, rather than Dirich- 
let conditions as in the standard case, and (ii) the resulting boundary layers are "weak" in 
the sense that they are O(e) amplitude, where e is the order of the viscosity. 

Although the existence of this new type of boundary layer, with quite different behavior 
from the standard type, is surprising to us, such layers have physical relevance (sse Appendix 
ICl) . In particular, one must understand these layers in order to treat cases arising in physical 
applications to suction-reduced drag in aerodynamics. Their analysis requires the study of 
hyperbolic boundary- value problems with Neumann or mixed Neumann-Dirichlet boundary 
conditions, an area that appears not to have received much attention, despite the extensive 
study of noncharacteristic hyperbolic boundary- value problems. We have described two 
approaches to these hyperbolic boundary problems, one involving a reduction to a problem 
with pseudodifferential Dirichlet conditions, and the other involving a reduction to a Cauchy 
problem on the boundary. We have provided examples where each approach works, but 
much work remains to be done on the general case. An important example where the second 
approach works is the case of supersonic inflow for the full compressible Euler equations 
considered in Appendix [Cl 

To study the small viscosity limit in the quasilinear hyperbolic-parabolic boundary prob- 
lems considered here, our approach requires estimates for the linearization of the problem 
about an approximate solution. The derivation of such estimates is completely open for 
cases other than the pure Neumann totally incoming case treated in the remainder of the 
paper. 

2 The quasilinear totally incoming case 

We turn now to our main task, the full treatment of the quasilinear case with full 
Neumann boundary conditions and totally incoming modes. 
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2.1 Construction of an approximate solution 



By plugging u® as in (|1.32p into the boundary problem (|1.3ip . Taylor expanding Aj(u*) 
about uo, and equating coefficients of equal powers of e on right and left, we obtain the 
following sequence of boundary problems: 



(°) ^2 A j( u o)djU = f, d d u \ Xd=0 = 

3=0 

d d 

(b) } j A j (u )djUi + } y dgAj (u ) {ui , dju ) = Au , d d ui\ Xd=0 = 

(2.1) 

d d y ' 

(c) ^ Aj (uo)djU 2 + ^ d u Aj (u )(u 2 , djup) = 

d d 

Am - ^ d u Aj(u )(ui, djUt) - y~] d 2 u Aj (n ) (ux , u x , djup) , d d ui\ Xd=0 = 
3=1 i=i 

and so on, where / = in t < and Uj = in t < for all j. Here / G -£P(R^ +1 ) for s 
large to be specified later. 

To solve (|2.ip (a). we first solve the symmetric, pure initial value problem on x d = 0: 

d-l 

Y J A M d i v = f\x d =o, w = 0int<0, (2.2) 

3=0 

and then the symmetric, dissipative boundary problem on £It f° r To small: 

d 

^A j (up)d j up = f, u \ Xd=0 = v, n = 0int<0. (2.3) 

3=0 

From (|2.2p . fj2.3[) . and the invertibility of A d {up) we obtain d d up\ Xd= o = 0. The subsequent 
linear problems (|2.ip (b).(c)....for the unknowns u\, u 2 ,-.- are solved by the same method. 
Standard theory gives < To < T\ such that|£| 

v G H s -\bn Tl ), n G H s -\n To ), u x G F s " 4 (fi To ), n 2 G ^ 7 (^t ), , „ 

m^-l-Sfc/o \ 2 - 4 

...,u k eH s 1 3k {n To ). 

Moreover, as long as s — 3M — 2 > 2±1 ; it i s eaS y t Q check that the remainder R £ in (jl.33p 
belongs to H s ~ 3M ~^(Q,t )- We now summarize this construction. 

Proposition 2.1 (Approximate solutions). Fix M G N. Consider the boundary problem 

(fL31~j) . u^ere / G £P(!+ +1 ) /or some s > 3M + 2 + Then (fOTj) /ias an approximate 

solution of the form 

u a £ (x) = u (x) + em{x) + ■ ■ ■ + e M u M (x), (2.5) 
satisfying (|1.33[) . u>/iere n^ G ^P~ 1_3fc (f2T ) an <^ remainder R e G -£P~ 3M ~ 3 (f2;r ). 

9 The drop by three units of regularity at each stage is due application of the Laplacian and the taking 
of a trace. Here we have chosen to restrict the Sobolev indices to lie in N. 
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2.2 Error equation 

We look for an exact solution of the form 

u £ = u a + e L v £ , where 1 < L < M. (2.6) 

To obtain the problem satisfied by v we divide the equation £{u) — £(u a ) = —e M R £ by e L 
to obtain 

d 

Aj(u a + e L v)d Xj v + E{u a , Vu a , e L v)v - eAv = -e M - L R £ 

3=0 

d Xd v\ Xd =o = 
v = in xq < 0, 



(2.7) 



where with V = (d Xl , . . . , d. 



E(u a ,Vu a ,e L v)v :=Y^(^f o d u Aj(u a + se L v) ■ v ds^j d Xj u a . (2.8) 
To obtain a linear operator acting on v on the left we rewrite (|2.7[) as 

d 

^ Aj(u a )d Xj v + E(u a , Vu a , 0)v - eAv = 

3=0 

- e M - L R e + e L B 1 {u a ,e L v){v, Vv) + e L B 2 {u a , Vu a ,e L v){v, v) := F e (v, Vv) ( 2 - 9 ) 
dx d v\ Xd=0 = 
v = in xo < 0. 

Here B\ and B 2 , defined by the equation, are smooth functions and bilinear in their last 
two arguments. 

Next we rewrite (|2.9p as a 2N x 2N first-order system for the unknown U = (ui, i^)* := 
(v,ed Xd v) 1 , setting d" = (d Xl , . . . , d Xd _ x ): 



where 



d Xd U = ~G(p(x),ed x ,)U + F £ (U,d"u) 

TU := u 2 = on Xd = 
U = in x < 0, 



F - {U) ( -Mv,Vv)) andG(p(x), E ^) = (^ with 



(2.10) 



d-i (2.11) 
M = ^ A j {u a )ed X] + sE(u a ,Vu a , 0) - e 2 A x » and 4 d = AK* 
3=0 



lu a ) 
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In (|2.10p we have set 

p(x) = {pi(x),p 2 {x),p 3 (x)) where 

pi(x):=u , p 2 (x) = u a -u , p 3 (x) = eE(u a ,Vu a ,0). 

To prove weighted estimates we introduce U = e~' yx °U, F = e~~ <X0 F, where 7 > 1, and 
observe that ()2.10f) is equivalent to 

d Xd U = -G^(p(x) 1 ed xl ,e 1 )U + F e (U, d"U), 

TU := u 2 = on x d = ( 2 - 13 ) 
U = in x < 0, 

where G 7 is defined by replacing d XQ by d XQ + 7 the definition of G. 
2.3 Symbolic preparation 

The operator G 7 in (12.13P is the semiclassical differential operator defined by the symbol 
where, with p = (pi,p 2 ,£>3), /? = (Pa, ■ ■ -,Pd-i,l') 

d-1 d-1 

M (p, p) : =ip + 1 ' + J2 Aj (pi + p 2 )i/3, + p 3 + Yl $1 

3=1 3=i (2 ' 15) 

A(p) := A d ( Pl +p 2 ). 

Lemma 2.2. For p\ £U, {p 2 ,p 3 ) in a small enough neighborhood oo 2 x w 3 o/(0, 0), and P 
in a small enough neighborhood up of 0, there exists a a C°° invertible matrix T(p, P) such 
that T~ 1 G(p, P)T has the block diagonal form 

T- 1 GT=(* ° p Y (2.16) 



where 



with 



and 



T <P>n=l-A-iM + n /+r 2 )' (2 - 17) 



r 1 (p,p) = (0(P) + 0( P3 )) 2 , r 2 (p,p)=0(P) + 0(p 3 ). (2.18) 



H(p,p) = -A~ X M + t\ 
P( P ,P) = A + At 2 . 



(2.19) 



Proof. The proof is a simple computation. Look for T of the given form and use the 
invertibility of A to solve for n, t 2 by contraction. 

□ 
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2.4 Computation of the low frequency Evans function 

Consider the N x N parabolic problem 

d 



d t u £ + ^2 Aj(u £ )djU £ - eAu £ = f 

1 (2.20) 

d x u e \ x = = 

u £ \ t< o = 



We now examine the Fourier-Laplace transform of the linearization of (|2.20p about a con- 
stant state u = u G U, where IA is the neighborhood of specified in Assumption 11.51 
Writing Q = (r, rj, 7) for now and setting Aj = Aj(u) and A(irf) := Yli~ l Ajirjj, we obtain: 

(it + i)v + A(ir])v + A d v Xd + e\r]\ 2 v - ev XdXd = f 
v Xd = on x d = 0. 

By multiplying through by e and rescaling x d and frequencies {x d — > ^f, ( — > e() we reduce 
to the case s = 1. Rewriting (|2.21|) as a first order system we obtain with U = (u 1 ,^ 2 )* := 



TJ7 = U2 = on Xd = 0, 



(2.22) 



where 



F ~ {-ef) ' G(0 : " (<r + 7 + Aiirj) + M 2 X) " ^ 

For £ 7^ let E~ be the stable generalized eigenspace of G(C)- Define the Evans 
function 

D(C) = det(£-(C),kerr). (2.24) 

Nonvanishing of the high frequency Evans function (a rescaled version of D (£)) was verified 
in IGMWZBj . Prop. 3.8. For fixed < r < R the fact that D(() ^ for r < \(\ < R is 
proved in section 4.1 of 

(GMWZ5]0rhus, 

we focus now on the low frequency region. 
We show that the Evans function vanishes in the limit as £ — > 0. For |£| small we 
conjugate G(C) to a block diagonal form 

S-HCMOSiC) = (^ C) p ° c) ) :=G H , P , (2.26) 



10 More precisely, the estimates (4.7) and (4.8) in IGMWZ5] are also true with 5RA replaced by |A| on the 
left. Those estimates and Sobolev's inequality readily imply the trace estimate 

|«(0)| < C{r,R)\v x (0)\ (2.25) 

for (v(0), v x (0)) 6 E~(Q and £ in this frequency range. 
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where 

H(() = -A" 1 (ir + 7 + A(ir,)) + 0(p 2 ) (p = |C|), P(C) = A d + O(p), (2.27) 
and the conjugator can be chosen to have the form 

5(C) = ( a * Sl l {C) ) with 5 21 (C) = O(p). (2.28) 



5 2 i (C) I 

To construct 5 one can simply look for a matrix of the form (|2.28[) satisfying GS = SGh,p, 
and use the invertibility of A d to solve for the off-diagonal blocks of 5 and the error terms 
in (12371) . 

Writing GS = SGh,p and equating (1,1) entries we obtain 

5 2 i (0 = H((). (2.29) 

Set U = S(()U, where U := ( UH ) and consider the equivalent problem 



(2.30) 



up 

U XA = G H ,pU + S^F 

f(C)U :=TS{QU = H{Qu H + u P . 



Let F"(C) = S-\QEr{C). Since A d is positive, F~(Q = {(z,0) : z G C^}. On the 
other hand we have from f|2.30|) 

kerf(C) = {(it;, -H(()w) : w G C^}. (2.31) 

This gives immediately 

D{() = det(F-(C),kerf(C)) = detff(C) for p small, (2.32) 
where each equality holds up to a factor that remains bounded away from zero for p small. 

2.5 Resolvent estimates by degenerate symmetrizers 

Recall that F~(() = {U = (u H ,0) : u H G C^}. Thus, for U G F~(C) we have 

f (C)W = H(()u H , (2.33) 

so 

\U\ = \u H \ = \H-\0HCM (2.34) 
This gives the degenerate trace estimate 

|f (C)W| > R(0M: wh ere R(() := \H-\()\-\ for U G F~{Q. (2.35) 
Proposition 2.3. Let p := Then for r > small enough we have 

\R(0\ >C(j + p 2 ) for 0< |C| <r. (2.36) 



llr This Proposition does not require Ad > 0; it remains true when H satisfies the generalized block structure 
property of [GMWZ6] . 
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Proof. 1. Write H{() = pH((, p) and fix C G S+. For (C, p) in a neighborhood of (£, 0) we 
use the smooth block reduction of ( [GMWZBj . (3.20)) 

V~ l HV = diag(^ fc ), (2.37) 

where H k has spectrum in a small disk centered at for p k the fcth distinct eigenvalue of 
H(C,0). By compactness of it suffices to show 

1^(01 < (2.38) 

for ( in a neighborhood of any fixed C £ st_ and p small. 

2. Let E be a Kreiss symmetrizer constructed as in [GMWZ6J for H k . The symmetrizer 
E satisfies 

(a)»(E^ fc )>C(7 + p), 

(6) |S| < C (2 ' 39j 

near the basepoint. The estimate (|2.39p (a) implies that Ei5& is invertible near the basepoint 
for 7 + p > 0, and since the same is true for H k , we see that E itelf is invertible near the 
basepoint for 7 + p > 0. The estimate (|2.39p (a) also implies 



C{l + p)\u\ 2 < $t(ZH k u,u) < \EH k u\\u\, (2.40) 



so 



KEFfc)- 1 ! < -SL- (2.41) 
7 + P 



and thus (since E is invertible) 

a 



H k 1 \ = |(EFfc) _1 E| < — — . (2.42) 

7 + P 

□ 



2.5.1 Resolvent estimates 

With U = as in (jXHD , we have 



^"(0 = {(«H,0) : « H € C N }, F+(C) = {(0,« P ) : up G C JV }, (2.43) 

where i ?=F (C) is the negative (resp. positive) generalized eigenspace of Gh,p(C)- Writing U 
as U now, we consider the problem 

d Xd U = G H , P U + F 

(2.44 

new = 9. 

Let |ujj|2 denote the L 2 [0,oo) norm, and let |it| be the norm of the trace at Xd = 0. 
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Proposition 2.4. Fix r > small. For < \Q < r we have the following estimate for 
solutions of (|2.44p .- 

(7 + P 2 ) 3 Wh\ 2 2 + \u P \l + ( 7 + P 2 ?\uh\ 2 + \u P \ 2 < C {\Fp\l + ( 7 + p 2 )\F H \ 2 2 + \g\ 2 ) . 

(2.45) 

Proof. 1. We use a degenerate symmetrizer of the form 

*(o-( ( ' + T' (0 

where /c > will be chosen sufficiently large, Sh(C) is a standard symmetrizer for the -ff(C) 
block (constructed as in [MZ1] , e.g.) and satisfies 

K(S H H) > (7(7 + /> 2 ) (2.47) 

while Sp satisfies 

Op — op 

M(SpP) > I (2.48) 
Spup ■ up > |up| 2 . 

Taking the real part of the L 2 [0,oo) inner product, (•,•), of —S^U with (|2.44p and 
integrating by parts gives 

~S fc l7(0) • C/(0) + (E7, U(S k G H ,p)U) = R(-S k U, F), (2.49) 

so 

\ (%p(0)| 2 - ( 7 + /> 2 )>h(0)| 2 ) + (7 + p 2 ?\u H \l + k\u P \ 2 < 

\(7 + p 2 ) 2 S H u H ,F H )\+k\(S P up,Fp)\ < ( 2 - 5 °) 

5( 7 + p 2 f\u H \l + C«( 7 + /9 2 )|F//|2 + <5%p|| + C^IFpI 2 . 

After absorbing interior terms in the obvious way from the right, it remains only to estimate 
the boundary terms. 

2. Using (|2.35p and (|2,36p . we have for the boundary terms, 

k\u P (0)\ 2 - (7 + P 2 ) 2 \uh(0)\ 2 = k\u P (0)\ 2 + (7 + p 2 ) 2 |up(0)| 2 - 2(7 + p 2 ) 2 \u H (0)\ 2 > 

2 

|2 1 / 1 „2\2i . /n\|2 



A ; |up(0)|^ + (7 + / 9^|uh(0)| / -C 



r(C) 



> 



Mo) 


%p(0)| 2 + (7 + P 2 ) 2 |^(0)| 2 - C\g\ 2 - C\u P (0)\ 2 . 

(2.51) 

For k large enough (12.51 j) and (|2.50|) imply the estimate (I2.45H . 

□ 
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2.6 The basic variable coefficient L 2 estimate 

Notation 2.5. 1. For u{x) G L 2 (R+ , iP (Rf, ) ) and ( = (C',7) = ( Co,C",7), set 

l«l«,7 = KC) s w(C / ,xrf)| i 2 {c , ia , d) . 

2. For u(x') G H s (R d ) set (u) s = \ {() s u\ L 2 (c) . 

3. Let A(eC) = (1 + (ej) 2 + (eCo) 2 + K"| 4 )^- For u ( x ), v(x') set 

\u\ A = \A(e()u((',x d )\ L 2 {( , tXd) , {v) A = |A«)w(C')Il2(c')' 

and similarly define |uL, (v)j, for other weights 4> = 4>(£,C)- 
4- For u(x) set (u)^ = (u(x' , 0))^. 

For given p(x), F, and g we now consider the following linear boundary problem corre- 
sponding to (|2,13p . where now we drop tildes and the superscript 7 on G: 

d Xd U--G{p(x),ed xl ,e 1 )U = F 

TU = g on x d = ( 2 - 52 ) 
U = in xo < : 

Our goal is to prove the following (degenerate) L? estimate for solutions of (|2,52p . 

Theorem 2.6 (Main 1? estimate). Under Assumption \1.5l there exist positive constants 
C , eo, 70 such that for all 7 > 70, < e < eo with < 1, solutions to (12.52P satisfy 

e\U\ Q + e(U) < C [y/i\F\ + (g) ) . (2.53) 

The preceding estimate is a composite of three more precise estimates corresponding to 
the three natural frequency regimes in the problem, the regimes in which e£ is of small, 
medium, or large size. 

Recall /3 = ((3',^') el^x R + is a placeholder for eC,. We shall localize with respect to 
the size of j3 using smooth cutoff functions Xj(/3), j = S, M, L, such that 

Xs(P) + Xm(P) + Xl(P) = 1, (2.54) 
where for some constants R\ (sufficiently small), R2 (sufficiently large) 

supp xs C {0 < |/3| < 
3 

supp xm C {-R\ < \P\ < R2] (2.55) 
3 

supp xl C {-R 2 < 

Notation 2.7. 1. We will occasionally use the symbol xm to denote a cutoff distinct from 
the one in (|2.55p . but also supported in a bounded region strictly away from the origin. 
Similar statements apply as well to xs> XL- 
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2. Choose smooth cutoffs Xi(P)> X2(P) identically equal to 1 near j3 = and compactly 
supported in ojr such that 

XiX2 = Xi, XsXi=XS- (2.56) 

3. The symbol ro will always denote a symbol or operator of order zero. 

4- Denote by O(eD) a semiclassical operator with symbol s(x,f3) such that s = (3- f(x,/3) 
for some smooth /o 0(e) denotes an operator with symbol s = ef(x,(3) £ Soa. 

In a similar way define 0(e 2 ), 0((eD) 2 ), etc.. When speaking of symbols instead of 
operators we'll use, as before, the notation 0(eQ), 0(e), etc.. In ambiguous cases like 0(e), 
the intent (symbol or operator) should be clear from the context. 

5. Write the solution to (|2,52[) as U = (u,v). Define 

U A = (Au,v), (2.57) 

where A(eD) is the multiplier associated to the symbol defined in Notation \2.5[ 

Here are the estimates by frequency size: 

Proposition 2.8. Using the notation just introduced, we have the following estimates for 
solutions to (12. 52ft . Let R± , R2 be as in (12. 55ft . For R\ sufficiently small and R2 sufficiently 
large, there exist constants C , 71, ei such that for all 7 > 71, < e < e\ with e~f < 1 

(a) \xs,DU\^3 + j p3 + (xs,dU) £1+£ 2 p 2 < 



C V^\F\ + (g) + e\U\ + \x2,dU\ a + \xm,dU\ + e(U) 

(b) \xm,dU\ + V^{xm,dU) < C (e\F\ + y/e{g) + e\U\ + e(U) ) 

(c) |Xi,-D^AlvA + Vc(xl,dU a )o < C (e|F| A -i /2 + Ve(g) + e\U A \ A -i/ 2 + e(U A ) A -i /2 ) ■ 

(2.58) 

Proof. The estimates (j2.8|> (b) .(c) are proved in |MZ1] . In the latter case we have applied 
the high frequency estimate of Proposition 4.6 of [MZl] after commuting (A~ 1 / 2 )d through 
the problem. We concentrate now on proving (j2.8|) (a). 

a. Localize to small frequency region. Commuting xs,D through (|2.52|) . we see 

that xs,dU satisfies 

Xs,DU Xd - -G d xs,dU = Xs,dF + -[xs,d,G d ]U 

e e (2.59) 

Fxs,dU = xs,D9 on x d = 0. 

There is a high frequency contribution to the commutator because of the x' dependence of 
G, and to get a good estimate for this we use the semiclassical calculusEl Since 

Xs,dG d = (xsG) D + -.(dp X sd x iG) D + e 2 r , (2.60) 



12 Since s must be bounded, we must then have |/| = 0(1/|/3|) for |/3| large. 

13 Even though the symbol of G is not bounded, one can use and directly estimate the formula for the 



remainder given in (A.6) of [GMWZ2] to prove ((2760]) . 
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we have 



\[XS,D,G D )U = ^(dp,xsd x iG) D U + er Q U. (2.61) 
Thus U a = xs,dU satisfies 

d Xd U a - -G D U a = F a (2 g2) 
FU a = g a on x d = 0, 

where 

\F a \o<C\F\ + \(dpxsro) D U\o + e\U\ , (g a } < (g} . (2.63) 

To prove (|2.58p (a) it suffices to prove the same estimate with xs,dU, F, and g replaced 
by U a , F a and g a . 

b. Conjugate to Ghp,d- Let T(p,f3) be the conjugator constructed in Lemma 1231 
and set 

°hp = ( H n p ° ) (2-64) 
\0 P + er J 

Extend T(p(x),(3) smoothly to all j3 G M. d x IR + as a semiclassical symbol with a uniformly 
bounded inverse, and use the calculus to construct right and left (approximate) inverses 
T_i 5 £) satisfying 

TDr -^ = / + t (2-65) 
T_i, D T D =I + e r . 

The right and left inverses are not equal, but we use the same notation for both. The 
symbol T_i in each case has the form 

T_i(p(x),/3) =T- 1 + er . (2.66) 

Defining V = T-i^Ua, we have 

(a) T D V = U a + e 2 r U a 

(b) (d Xd T D )V + T D d Xd V = d Xd U a + O(e)(r U a + eF a ) = (2 _ g7) 

jG D T D V + F a + O(e)(r U a + eF a ). 
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We have the following symbol equalities 

(") T =(i A ^+O(e0 + O(e) 
(6) T-i = (I + 0(e() + 0(e) 

(c) T^d Xd T = (° r Q °) + O(eC) + 0(e) 

(d) GT X2 = (° 2) X2 (eC) + O(eC) + 0(e) ^ 

(e) edpT.x = eO(eC) + O(e) 

(/) l -(ed p ,T- l )d x ,(GT X 2) = (° ^ X2«) + 0(e() + O(e) 
( 5 ) -T_!GT X2 = -G HPX 2 + f? X2«) + O(eC) + 0(e). 



e e V° r o. 

For dZ68|)(g) we used 03322), (EHI, and (l238]l fd N ). 

Applying the operator T_i £> to (|2.67p (b) and using the semiclassical calculus, we obtain 
in view of the symbol equalities ([2.680 : 

^ = \ { H 0° PD + er ) X2 > dV + roFa + ° {e)Ua + ° {eD)V + ° {e)V - (2 - 69) 

Observe that terms on the right in (|2.68j) (c).(f). and (g) all make contributions to the ro 
entries of the first matrix on the right in (|2.69j) . Using the calculus to commute Xl,D through 
(|2.69p . we obtain 

d Xd (xi, D V) = \ ^ p D 7 ero ) (xi,dV) + r F a + 0(e)U + 0(eD) X i,dU + (r fyxi) D l7. 

(2.70) 

Next define 

F b := r F a + 0(e)U + 0(eD) X i, D U + (t^x^dU, (2.71) 
and observe that since U a = Te,V — e 2 roU a = TdXi,dV + s 2 roU and Ub = Xi,dV satisfies 

d Xd U b = ±( H ° p er 2 )u b + F b 

e \ P D + er J (2.72) 

TT D U b = g a + e 2 r U := g b on x d = 0, 

to prove (|2.58p (a) it now suffices to prove the same estimate with xs,dU, F, and g replaced 
by U b , F b and g b . Observe that \feF b is a sum of terms including ^0(eD)xi,dU . The 
latter term is absorbed using the following Lemma, whose proof is elementary. 
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Lemma 2.9. Fix 5 > 0. Then for 7 large we have 

\ 3 / 1 5 O 

(1) £2^ < (J I £^2 +£2p J 

(2) eV < <5 

Define 



(2.73) 



p £ ;°,. o ). (2.74) 

A direct computation using the invertibility of P shows that for (3 £ ujp one can choose a 
matrix symbol T c of the form 

T e (p{*),0) = ( I o 7) (2 ' 75) 



T- l G b T c ={ H n _ ° ) =G HP . (2.76) 



such that 



P + er 

As before we extend and invert T c ,£>- The operator T c ,-i £) associated to the symbol 

^-i=(j 7°) (2.77) 

is easily seen to be a right and left inverse satisfying the analogue of (|2,65D . 

Redefine V = T c _\^dU\,. Now repeat the preceding argument line for line, but note, for 
example, that instead of (|2.68|) (c).(e).(f) we have, respectively, 

T c ^d Xd T c = 0(e) 

^ ,Tc '- 1= (o C 0°) (2-78) 
^(ed l3 >T c ^ 1 )d x/ (G b T c ) = 0(e). 

We set U c = Xi,dV and use the calculus just as before to find that U c satisfies 

1 

~ e L ^'— ' (2.79) 
(b) TT D T CtD U c = g a + e 2 r U := g c on x d = 0, 

where F c has a formula like (|2.7ip (with P5 in place of F a ). Thus, to prove (|2.58|) (a) it now 
suffices to prove the same estimate with xs,dU , F, and g replaced by U c , F c and g c . 

c. Block structure. Recall that Ghp is given by (I2.64p . where H(p,f3) and P(p,/3) 
are as in LemmaO Let p' = (pi,p 2 ), define H(j/,P) = H(rf,0,p), V{p',j3) = P(p',0 5( 3), 
and set 



(a) a^C/c = -G H p,dU c + P c 
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Note that for f3 G ojp 

(H(p,/3) \ (r Jn,(0(pa) \ 

I P(p,)9)J- fop(p ^ )+ i 0(p 3 )J' (2 ' 81) 

and thus 

G HP (p(x),P) = Ghp(p'(x),P) + ( £r Q ° £ °J . (2.82) 

To proceed further we need to conjugate Qhp to block structure form, which is especially 
simple in the totally incoming case. Introduce polar coordinates 

P = f/j3, where /3 G 5"^ = {(/?', 9) G S d : 9 > 0}, p' = |/3| (2.83) 

and write 

n(p',p)=p"H(p'J,p'). (2.84) 

Similarly we set ( = (C>7) = C/ICI an d p = Id- 
Proposition 2.10 (Block structure). Let p' £U. For each (3 G there is a neighborhood 

O of (?/,/?, 0) in R 2N xSf xl + and a C°° matrix V(p',/3,p') defined on O such that 
V^ 1 T~LV has the following block diagonal structure: 

1. 7/9 > 0, then V- X UV = Q where = (Q + Q*)/2 < c < 0. 

2. When 7' = 0, we have 



v- l nv 



(p'J,p'):=h(p'J,p'), (2.85) 



where the qj are scalars, not necessarily distinct, such that $tqj = when 7' = p' = 0, 
d~,(9tqj) < c < and d p i{^Rqj) < c < 0. 

There is a C°° matrix W(p',f3) defined on a neighborhood of (p',0) such that 

W~ l VW = p(p',/3), where Kp > C v > 0. (2.86) 

Proof. A general block structure result that applies in our case is Lemma 2.10 of [MZlj . 
The simplification due to the totally incoming assumption, A2 > 0, is explained in Corollary 
7.9 of (GMWZ6]0 □ 

d. Degenerate symmetrizers. The simple block structure described in Proposition 
12.101 permits the following simple construction of degenerate symmetrizers. Let f2-r = {x G 
l+ +d : < x < T} and bQ T = {x G tt T : x d = 0}. 



14 This assumption rules out glancing modes, and also guarantees that all blocks are lxl near points 
where 7' = 0. 
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Proposition 2.11. Fix x G bQ? o-nd ( G S 1 ^ and consider a neighborhood O of (p'(x), 0) 
£ + on which a conjugator V(p' , (3, p') as in Proposition ^. 1(X is defined. For 



S h 

Sr 



in R 2N x Sf x 

(x,C) suc/i i/iat (p' (x) , ( , e p) G 0, define 

S(C,ep) 

where the N x N matrices Sh, S p are given by 

S h = -{e 2 1 2 + e A p A )I N , S P = KI N , K > 0, 
and set h(p' , (3 , p') := p'h.(p' , (3, p'). Then, depending on O, either 

$t-Sh{(,£p)Hp'{ x ),(, £ P) = P{ £2 1 2 + £ 4 p 4 )k{x,(,ep), where k > C > 



or 



£Y+£V)(7io,i+£A,i: 



(eV + ^VXT&cyv + epXiv). 



(2.87) 



(2.88) 



(2.89) 



(2.90) 



where boj(x,(, £ p) > C > 0, bi t j(x,£,ep) > C > 0. Ako, 

sft-5 p p = -K %tp(j/(x),0) > ^KC P , where K > 0, C p > 0. 

Finally, for u = (uh,u p ) G C 2JV we have 

(S((,£ P )u,u) > K\u p \ 2 - ( £ 2 7 2 + £V)K| 2 - 



(2.91) 



(2.92) 



Proof. The equalities (|2.9ip and (12.92H are immediate. Consider O as in case 2 of Proposi- 
tion 12.101 The properties of qj stated there imply 



qj{p', (3, p) = 7'aoj + pa ltj + idj 



(2.93) 



where ao,j, aij, dj are real functions of (p',/3,p') such that aoj < c < 0, aij < c < 0. 
Setting f3 = eC, and noting that (3 = ( and p' = ep, we see that (|2.90p holds with 

b ,j(x,(,ep) = -ao,j{p'{x),(,e(), bij(x,(,ep) = -ai tj (p(x), (, e(). (2.94) 

Similarly, (|2.89|) holds when O is as in case 1 of Proposition I2.1U1 

□ 

e. Microlocalize. Next we construct a pseudodifferential partition of unity that will 
allow us to prove estimates using the symmetrizers just constructed. 

Let K p i be the compact set given by the closure of the range of p'(x) on ^tEI For 5 > 
small we can choose an open cover {Ok} of K p i x Si x [0, S] by sets Ok on which conjugators as 



15 Recall that, by finite propagation speed of the hyperbolic problem, p(x) is constant outside a compact 
subset of Qt- 
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in Proposition 12.101 are denned. Choose partitions of unity m(x) and i/j m (j3,p') subordinate 
to open covers of Qt and Si x [0, fi\, respectively, with the property that for any pair (I, m) 
there exists a k such that 

(x,/3,//) e supp^(x)V m (/3,p') {p'(x)J,p') G O fc . (2.95) 
Next define the bounded families of classical symbols 

#,mO>0 = Ki(x)ip m (C,ep)- (2-96) 
After re-indexing this family as <^(x,£), we rewrite the unknown U c in (|2.79|) as the finite 



sum 



U C = J2 U h wh ere U t := </>f )0 77 c . (2.97) 
i 

Next we commute $\ y) through the problem (|2.79p . Observe that (|2.79p (a) is unchanged 
if Ghp,d is replaced by Ghp,d X2,d, and henceforth we include the factor \2 (often sup- 
pressed) in the definitions of H and V (|2.80p . Observe that Vx2 is a classical symbol of 
order zero and that 

H = epn(p'(x),(,£ P )x2(e() = eU\ (2.98) 

where H* is a classical symbol of order one. The leading terms in the symbols of the 
commutators [H* D , (pf D ] and -\Pd,</>i £>] are, respectively, 

(a) d ( ,U* D x «\)\ - 9 f /0f D X >H* e C° 

1 1 (2 99) 

(6) - (0 f ,p D x ,<j>! - d c ,<j>i D X ,V) G - C" 1 . 

Thus, from d2321) and ((275)1 (a) we find: 

Scl/j = ^hp,dC/z + (g i r °_ J C/ c + 0f iD F c (2.100) 

where r_i G C _1 . 

The boundary operator TTf)T Cj £) can be viewed as an element of C°. Its leading symbol 

is 

4/ + „ /+ r 2 ) (o ?) = (« + £r » ' + £r » + «») ' 

(2.101) 

Using (|2.99p (a) (and a similar computation for the commutator [0(eD)x2,Di <t>f dY) we ^nd 

TT D T c , D Ui = (f>l D g c + er U c . (2.102) 

For a fixed Z, set 

U c = (U h , U p f, Ut = <f>l D U c = (u h , up)*, and tf jD F e = (f h , f p ) 1 . (2.103) 



TTT, = (0 I) , ^ 
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We can rewrite (|2.100p and (|2.102p as follows: 

(a) c^u/, = -UdUh + r U h + f h 

{b)d Xd u p = -V D u p + r -^-U p + f p , r^GCT 1 ( 2 ' 104 ) 

(c) U D Uh + Id u p = r o9c + er U c . 

Here we have set I £ D = I + 0(eD)x2,d and used (|2.10ip and the semiclassical calculus to 
compute VTy)T c ^. 

f. Conjugate with Vd and Wd- For the same fixed I as in (|2.103j) suppose that for 
(x, Cj € supp <^f (x, C)> (p'C^)) C) e /°) i s contained in an open set O as in case 2 of Proposition 
I2.1QH . and let V(p'(x),(,ep) G C° be the corresponding conjugator as in (|2.85p . Extend 
V and (approximately) invert Vd in the classical calculus to obtain left and right inverses 
such that 

VdVS 1 = 1 + r_i, V^Vb = i" + r_i, where r_i G C -1 . (2.105) 
Defining = V L) 1 Uf l we obtain by a computation similar to (|2.67p (b) 

(d Xd V D )w h + V D d Xd w h = -UdVdWh + r £4 + r /ft. (2.106) 
Here we have used the fact that 

-Hd is a bounded family in C 1 . (2.107) 

Applying V^ 1 to (I2.1()(-jp and using (12. 107)) again, we find 

1 1 

d Xd Wh = -V^ 1 n D V D w h + r U h + r fh = ~^DW h + r U h + r f h , (2.108) 

where h = h(p'(x),£,ep) is as in (gSlD , and hence ~h.£> 6 C 1 . 

Similarly, extend the conjugator W(p',e/3) in (|2.86p and construct approximate inverses 
of Wd in the semiclassical calculus such that 

W D W D l = I + er , Wp l W D = I + er . (2.109) 

Defining w p = W^u p we obtain by using the semiclassical calculus and computing as 
above 17 ! 

d Xd w p = ^pdw p + ~^-U p + r f p , where r_i G C _1 . (2.110) 

g. Interior estimates. We quantize S/j and S p as in (|2.88p by setting 

S h , D = ~ (e 2 7 2 In + e\p 4 I N ) D ) X 2,d, S PtD = KI nX 2,d, (2.111) 



16 We omit the details for case 1, which is similar but easier. 

17 Here we have used roX2(eC) € -C _1 , which holds since £X2(eC) £ C _1 . 
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where (p a In)d £ C 4 . Pairing (|2. 1Q8[) and ()2.110p with Sh,DWh and S p> dw p , we obtain the 
identities 



(Sh,DW h ,w h ) + ^(S h)D h D w h , w h ) = -2U(r U h + r fh, S h)D w h ) 

1 r_i 
{S p ,DWp, Wp) + $t~(S p>D Pdw p , Wp) = -2K(— Up + r f p , S PiD w p ). 

Since 

\(rofh,Sh,DW h )\ = | (r A, (e 2 7 2 + £ 4 P 4 )d X2,DW h ) \ 

3 \ 

< 



(2.112) 



^i( £ 2 7 2 + D ro/fc> 1 (£ 2 7 2 + ^4)! X2 DWh 



(2.113) 



C5 e|Alo + ^l W fe| £7 3/2 +£ 5/2 p 3> 

we see that 

|2 , /-» „| r |2 , r|„„ |2 



-2»(r ?7 fc + ro/h,5 ft ,£,«;h) < C\U h \ 2 £1+£2p2 + C s B\f h \' + 5\w h \^ 3/2+e5/2p3 . (2.114) 



Similarly, since \r_iUp\ < %\U P \ , 



-2${(-±Up + r /p, S p ,dw p ) < — \U P \ 2 + C s e\fp\l + -\w p \ 2 . (2.115) 



Next set h* = \h 6 C 1 and note that 



KS^h^.^-CCSfth*)^,^)! <CeVKIo + Ce Kip"- ( 2 - 116 ) 
Here we have used the classical calculus to obtain, for example, 

(£Vx2«)) D h^ = (e 4 p\ 2 (eOh*) D + e 4 r 4jfl , r 4 e C 4 . (2.117) 
The semiclassical calculus implies X2,dPd = (X2P)d + stq, so 



~{Sp,D PDWp,W p ) - ^({S p p) D Wp,Wp) 



< C\w p \l. (2.118) 



Now we can use ()2.90p and the Garding inequality for the classical calculus to get 
estimates from below: 

^((S h h) D w h , w h ) > C{e 2 ^\w h \l + e 5 \w h \ 2 p3 ) - C{e 2 ^\U h \ 2 _ x + e 5 \U h \ 2 p2 ) u 

-C{eW\U h \l + e A l\U h \ 2 p ). 

To obtain (|2.119p we have used, for example, the Garding estimate: 

K {(e 5 p 6 b hjX 2(eC))DW h j,w h j) > Ce 5 \w h;j \ 2 p3 ) - Ce 5 \U h \ 2 p2 , (2.120) 

where b±j is as in (|2,90p and Wh,i 1S the j-th component of w^. The error terms in the 
second line of (|2.119p come from "cross-term" estimates like 

K ((e A 7P%,jX2(eO)DW htj ,w htj ) > Ce A 7 \w h>j \ 2 p2 ) - Ce^\U h \ 2 p , (2.121) 
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(2.123) 



Another application of the classical Garding inequality gives 

^((S p p) D w p ,w p ) > ^\w p \ 2 - j\U p \ 2 _ v (2.122) 
Combining the above estimates for Wh we obtain 

(eVKIo + £ 5 Klp3) + (S h)D w h ,w h ) < C s e\f h \l + C\U h \ 2 ei+e2p2 + 
C^\U h \\ + e 5 \U h \ 2 p2 + e 3 7 2 \U h \ 2 + e 4 7 |t4|*), 

after absorbing Wh norms from the right using Lemma 12.91 Similarly, we find 
K CK C 

— \w p \ 2 + {S PtD w p ,w p ) < C 5 e\f p \l + —\U p \l + -{UpW, (2.124) 

after absorbing w p norms from the right, 
h. Boundary terms. We clearly have 

{S h>D w h , w h ) > -C 2 {e 2 j 2 {w h ) 2 Q + e\w h ) 2 2 ), (2.125) 

and an application of the classical Garding inequality gives 

{S PjD w p ,w p ) > K(w p ) 2 - C(U P ) 2 _ V (2.126) 

We use the classical calculus and the fact that Hd E eC 1 to rewrite (|2.104p (c) as 

V^UdVdvjh + V^r D W D w p = r g c + er U c , (2.127) 

which implies with a new eroU c 

B D \ h ) := h D w h + l e D w p = r g c + er U c , where X £ D = V^I £ D W D . (2.128) 
\ w pJ 

Clearly, 

(BbB D h) , (Z h )) < C(9c)l + Ce 2 (U c ) 2 , (2.129) 



and we now proceed to estimate (B* d Bd ( h ) , ( h ) ) from below. 

\ w pJ \ w pJ 



Lemma 2.12. 



C^iwh) 2 + C x e A {w h ) 2 - C 3 (w p ) 2 - C {e 2 j 2 (U h ) 2 1 + e 4 (U h ) 2 + £ 3 7 <t4)o) • 

(2.130) 
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Proof. We start from 

(B* D B D (^) , (™*) ) = {h* D h D w h , w h ) + 2Vt(h D w h ,l £ D w p ) + \lhw p \ 2 > 

\ w pJ \ w pJ (2.131) 

- (h* D h D w h ,w h ) - C 3 (w p )l 

Using (|2.93|) we see that 

h = Hp'(x),(, ep) = diag(h j ) = diag(A, + iDj), (2.132) 
where Aj = ejaoj + e 2 p 2 a\j and Dj = epdj. Now 

h* }D h j;D = A* D Aj t D + D* }D D j;D + i{Al D D hD - D* >D A j)D ), (2.133) 

where 

i(Al D D j>D - D* D Aj, D ) = eVo + e 3 r 2 , r 2 E C 2 , (2.134) 
since, for example, the classical calculus implies 

(e 2 p 2 axj)* D (epdj) D - {epdj)* D {e 2 p % a Xd ) D = e 3 r 2 . (2.135) 

Next we compute 

A* jiD A jtD = (e 2 j 2 {a 2 0J ) D +e 2 j 2 r^ 1 ) + {e^clx^d + £ 4 r 3 ) + (2e 3 ~f{p 2 a 0J a hj ) D + e 3 jn) . 

(2.136) 

The classical Garding inequality and (12. 136f> imply 

( A *j,D A j,DW hJ ,w hJ ) > Cie 2 j 2 (w hJ )l + Cie 4 (w h j)l + Ce 3 j(w hJ )l 

- C (eV(^) 2 i + eHU h ) 2 + s^(U h ) 2 ) (2 . 137) 
- C (e 2 l 2 (w Kj ) 2 _ 1/2 + e A {w hJ ) 2 3/2 + e 3 -f{w hJ ) 2 1/2 ^ , 

where the error terms in the second line of ()2.137j) are Garding errors, while those in the 
third line arise from the composition errors in (|2.136|) . From (|2.133|) . (|2.134p . and (|2.137p 
we obtain, after absorbing error terms involving Wh and w p by taking 7 large: 

(h* D h D w h , w h ) > C l£ 2 7 2 K)o + C^iwh) 2 - C (eV(^)-i + e\U h ) 2 + eS(U h ) 2 ) . 

(2.138) 

Using (|2.13ip . and f|2. 138|) we obtain the estimate of the Lemma with new constants. □ 

Combining the estimates of this paragraph we find, for constants as in (12.125p . (I2.126p . 
Lemma |2. 121 and some M > to be chosen: 

{Sh,DW h , w h ) + {S p , D w p , w p ) + M(B* D B D ( Wh ) , ( Wh ) ) > 

\ w p/ \ w pJ 

(Md - C 2 ) (eVK>o + e\w h ) 2 2 ) + (K - MC 3 )(w p ) 2 {2A39) 
- MC (eV<£4> 2 -i + e\U h ) 2 + e 3 7 (^)g) - C{U P )\. 



40 



i. Conclusion. To finish the proof of Proposition 12. 81 we first add estimates (|2. 123j) and 
(|2.124p . and then add M {B* D Br>(whi w pfi (whiWpY) to both sides of the resulting inequality. 
Boundary terms on the left in the estimate so obtained are estimated from below using 
(12.1391) : on the right one uses (J2.129I) . After choosing M so that MC X > C 2 and then K 
such that K > MC3, we get (with a new C) 

(eVKIo + £ 5 Klp3 + jKlo) + (eV(™A>o + £ A {w h ) 2 2 + H)o) < 
C{e\f h \l + e\f p \l + (g c )l) + 

U h \% +e ^ + eWWnt, + e*\U h \l 2 + e^\U h \l + eHU h \ 2 p + ^-\V p \l + + 
C {eW{U h )\ + e\U h )\ + e^{U h )l + (U^ + e 2 (U c ) 2 ) , 

(2.140) 



where the last two lines are "error" terms. Since 

Uh = V D w h + r_iu h and u p = W D w p + er u p (2.141) 

for V D , W D as in (I2.105j) . (12.1091) . the estimate (12.1401) holds with {wh^Wp) 1 replaced by 
(uhjUpY = Ui. Recalling that 

U c = (U h , Upf = J2 U i and ^IdFc = Uh, fpf, (2.142) 



summing the estimates over I, and absorbing error terms from the right using Lemma |2.9[ 
we conclude 



eW\U h \l + e 5 \U h \l 3 + ^\U P \£) + (e^(U h )l + e\U h )l + (U p )%) < 



(2.143) 



C{e\F h \l + e\F p \l + {g, 



c/o) 



This estimate is stronger than the estimate described at the end of paragraph b as being 
sufficient to prove (|2.58p (a). This concludes the proof of Proposition 12.81 

□ 

2.7 Higher derivative estimates 

In this section we'll use the notation for norms introduced in section 12.61 We use d 
to denote some tangential derivative, one of do, . . . , dd-i- Sometimes dU will denote the 
tangential gradient of U, instead of just a single partial derivative of U. 

Notation 2.13. 1. For k = 1,2,... let U*> k = {{^) k U, {^f^dU, d k U). Here &U 
represents all possible tangential derivatives of U order j. 

2. Define U^ k simply by replacing U by Ua in the definition ofU* ,k . 
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Proposition 2.14. Under the assumptions of section 2, there exist positive constants C, 
€q, 70 such that for all 7 > 70, < e < eo with ej < 1, solutions to f|2.52|) satisfy 

\U*>% + (^*' fc )o < C + i^)o^ . (2.144) 

This follows immediately from the following more precise estimates. 

Proposition 2.15. Using the notation just introduced, we have the following estimates for 
solutions to (12.52j) , Let Ri, R2 be as in (12.55j) . For R\ sufficiently small and R2 sufficiently 
large, there exist constants C , ji, e\ such that for all 7 > 71, < e < e\ with ej < 1 



_ ^ 1 /,, TT*,k\ 

e-y+£ 2 p 2 — 



C (^|F*< fc | + (g*' k } + e|C/*' fe |o + \ X 2. D U*> k \ 3 2 2 + |xAf,D^*'*|o + ^* 

(6) IXM,D^*' fe |o + ^(XAf,Df/*' fe )o < C (e\F*'% + vW' fc >o + e\U*>% + e(U*' k ) 
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(c) \xl,dU**\^ + V~e(XL, D U*/)o < C (e|F*< fc | A - 1/2 + ^9*'% + e\U*/\ A - 1/2 + e<^ A; ) A _ 1/ , 

(2.145) 

Proof. The estimates in (b) and (c) follow directly from the higher derivative estimates of 
[MZlj in the medium and large freqency regions. These are estimates with 7 weights for 
the linearized problem, so one can simply apply them to the problems satisfied by for 
various j. 

As usual, therefore, we focus on the small frequency region. If we simply differentiate 
the equation and throw commutators on the right as forcing, those new forcing terms are 
too large to absorb in a straightforward way. To get around this problem we reprove L 2 
estimates for an appropriate enlarged system. 

1. Enlarging the system. We begin with a solution U of the linear system (I2.52p 

d d U - -GU = F 

YU = g on x d = ( 2 - 146 ) 
U = in x < : 

Let d denote one of do, . . . , <9d-i- Observe that (^U, dU) satisfies the enlarged system 

1. 



0, i €) - i m (t) - (§;) + 1 

on x d = 0, (2.147) 



dUj e\0 G)\dU)-\dFj^\z[d,G]{%U))> 



Y \dU \dg 



dU ) 



in xq < 0. 
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2. Localize to small frequency region. Let Xs( e C) be a small frequency cutoff as 
before. Commuting xs,D through ([2. 147[) we obtain (writing xs for xs,d) 



ddixsU*' 1 ) - \ °) (xsU*' 1 ) 



XS, 



G 
G 



(2.148) 



where 

|F'| < CdF*' 1 ^ + KdpXsW*' 1 ^ + ep^lo). (2.149) 

The second commutator was computed like the corresponding term in the previous section 
CM). 

The boundary condition is 

r o \ i _ i * 



() r , XsU*^ = XS5 1 '* • (2.150) 

The problem (|2.148p . (j2.150p can be treated just like (|2.62p . We may now repeat the 
argument of the previous section to obtain the desired estimate of U*> . Iteration completes 
the proof. 

□ 

f~U\ 7—1. 

Remark 2.16. // U*' 1 had been defined instead as ( J , the first commutator in (|2. 148|) 
would have produced an unacceptable Odf/*' 1 ^) error. 



2.8 Nonlinear stability 



Notation 2.17. 1. Recall |«|fe i7 = \{Cl u {Cj x d)\o- For k E N we have the equivalence of 
norms 

H*,7~ E l k ~ H \d a nW (2.151) 
|a|<fe 

2. Set \u\* = \u\l°° ■ 

3. Define 

IMU,7 = |«|fe,7 + |e«9u|A;, 7 . (2.152) 

^. Lei M and L < M be the positive integers appearing in the nonlinear error equa- 
tion (|2,9p . They can be taken arbitrarily large as long as the approximate solution u a is 
constructed with sufficiently many terms. 

5. 4>{^) always denotes an increasing function of '7. It may change from term to term. 

6. Set 9' = <&,..., d d -{). 
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We return to the nonlinear error equation (|2.13p . and again drop tildes and the super- 
script 7. Let k(xq) be a smooth cutoff which is identically one on [0,To]. We will solve 
(|2,13p on [0, To] using the following iteration scheme: 



8 d Un+l - -GU n+l = K (x )F e (U n ,d"U n ), 

TU n+ i = on x d = 0, 
U n+ i = in x < 0, 

where from (|2.9p and (|2.1ip we see that F e (U n , d"U n ) has the form 

F e (U n , d"U n ) = e L - 3 h(u a , Vu°, e L U n , e, e^°)(eU n , eU n 
+ e L - 3 f 2 (u a ,Vu a ,e L U n ,e,e'y xo )(eU n ,ed"U n ) 



(2.153) 



+ e M ' L R 



(2.154) 



s 

A + B + C, 



for smooth functions fx, fa. For ¥(U,d"U) := k(xq)F 6 (U, d"U) consider the nonlinear error 
equation 

d d U --GU = ¥(U,d"U), 
e 

TU = 0onx d = 0, (2.155) 
U = in x < 0. 

Theorem 2.18. Recall d is the number of space dimensions. Fix constants k,L,M satis- 
fying 

k-Z> d - 
2 

M-L- 2/c — i > 1 (2.156) 
L-3-2k-->l. 

Suppose the forcing term f in (jl.3ip is chosen in H s (M. d ^ 1 ), where s > 3M + 3 + A;, so that 
u a as constructed in Proposition \2.1\ yields a remainder R £ G -ff fc (^r )El Then there exist 
constants eq, 70 such that for all < e < eo, 7 > 70 satisfying 67 < 1, the error equation 
(|2.155p has a unique solution U satisfying the estimates 

||C/|U, 7 < e M " L - 2fc -^(7) 

\UU < 1 (2-157) 

\dU\* < 1 

for some 4>("f), an increasing function 0/7. 



s This is the same R E that appears in (12. 154|) . 
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Proof. The first few points are some preliminaries. 

1. Sobolev inequalities. 

For k — 3 > | we have 

(a) e\dUU < C( 7 )(e\U\ k . 2n + e\d d U\ k . 2a ) 

(b) e\UU<C( 7 )(e\U\ k . 3tl + e\d d U\ k ^ 7 ). 

2. Moser inequalities. 

For k S N let a = (a±, . . . ,a r ) with \a\ = a\ + • • • + a r < k, en G N. Suppose 
|^i|fc,7 + < oo. Then 
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\{d^ Vl ) ■ ■ • (a Q ^ r )| < cJ^hMll \viU) 



i=l j^i 

3. Relations between norms. Directly from the definitions we see 

{a)\U\ Kl < C\U*>% 

C (2-159) 
(b)\U*>% < ^\u\ Kl . 

Let Xl(cC) be a high frequency cutoff like the one in (12.145p (c). Observe that 

\\U\\ ka ~ |C/| fej7 + |xL(e0C/-)k 7 . (2-160) 

4. High frequency estimate. Here we make use of a slightly modified form of the 
high frequency estimate in (|2.145p (c) with g = 0: 

\XlU*/\ A + VeiXLU*/)^ < C (e\F*>% + e\U*/\ + e(U*/)o) . (2.161) 

We can absorb the high frequency pieces of U^ k in the two terms on the right in (|2.15|) (c) 
to obtain 

\xlU*/\ a < C (e\F*> k \ + e\U*> k \ + e(U*> k } ) , (2.162) 



and then use the main L 2 estimate (|2. 144j) to replace the right side of the above inequality 
by Cy/e\F*' k \ . When |eC| is large, we have ^ > C((). Thus, with d2T59|)(a) we may 
conclude 

\XL(edU)\ k „ < Cy/i\E*>\. (2.163) 

5. Induction assumption. Let the first iterate U\ be 0. Assume there exist £1(7), 71 
such that for < e < ei, 7 > 7x, and some ^(7) 

||^„IU, 7 <2e M - L - 2fc -^(7) 

\U n l < 1 (2.164) 

\dU n U < 1 
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The main step is to show, after decreasing e\ if necessary, that U n +i satisfies the same 
estimates. 

6. Estimate F n := ¥(U n ,d"U n ). Set A = k(x )A for A as in (12.1541) . and define B 
and C similarly. 

Applying the Moser inequalities we have 

|A| fc , 7 < C( 7 )e L - 2 |[4|fc, 7 > (2.165) 

where C( 7 ) depends on L°° norms of (u a ,Vu a ) and eU n . 

Write edU n = (1 — XL){tdU n ) + XL{cdU n ), and corresponding to this decomposition set 
B = Bi + B2. Since \e(\ < C on supp (1 — Xi( e C))> we have just as above 

|BiU i7 < C(j)e L - 2 \U n \k !T (2.166) 

For B2 we have 

|B 2 | fel7 < C( 7 )( e L - 2 |C/ n | fci7 + e L - z \ XL {edU n )\ Kl ). (2.167) 

Moreover, we have 

|C| fe)7 < m^ L - (2-168) 

Summing these estimates we obtain 

|F n | fe)7 < C( 7 )(e L - 2 |CZ„| fci7 + e L - 3 \x L (edU n )\ k ^) + e M - L </>( 7 ). (2.169) 

7. Estimate \\U n +i ||fc, 7 - In view of the main estimate ()2.144p . (|2.159p . and (|2.169p 

we have 

C C 
|0»+iU, 7 < C|^+ilo < -^R fe |o < ^r|F»U, 7 

From (12.1631) and (|2.169D we obtain 



\ X L{edU n+l )\ k]1 < C\¥*^\ < ^|F n | fci 

< C( 7 )(e L - 2 - 2k \U n \ k ,, + z L -*- 2k \XL(edU n )\ kn ) + e M ~ L - 2k 4>{l). 
Adding the previous two estimates we find 

||CWi[|fc, 7 < e L - 3 - 2fc -^( 7 )||C/ n || fci7 + e M - L " 2fc -^( 7 ). (2.172) 

Provided ei( 7 ) is chosen so that e L ~ 3 ~ 2fc ~2 C( 7 ) < |, the induction assumption and 



(I2T72D imply 



\U n+1 \\ ka < 2e M - L - 2k ~hh)- (2-173) 
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8. L°° estimates. The equation gives 

£\ddU n+ i\k-2,j < C\U n+ i\k n + e|F n |fc_ 2)7 . (2.174) 

From ([2.169P we get 

|F„| fe , 7 < e L - 3 C( 7 )\\U n \\ k ^ + e M - i 0( 7 ). (2.175) 



Thus, 



e\d d U n+1 \ k ^ < 2e A// - i - 2fc -^( 7 ). (2.176) 



This together with the inequalities (|2.158p and the assumption (|2.156p immediately implies 
that for ei small enough 

f»n ll f< < 2 - 177 > 

e\dU n+1 \* < e. 

This completes the inductive step. 

9. Contraction. Thus, the sequence of iterates satisfies the estimates (12.1640 . One 
can now consider the problem satisfied by U n +i — U n and use estimates like those above 
(but simpler) to show that for ei small enough, the sequence converges to some U in the 
1 1 0,7 norm. A standard argument (involving interpolation and weak convergence) implies 
that U solves the error equation (I2.155P and satisfies the estimates (|2.157p in Theorem 12. 181 

This completes the proof of Theorem 12.181 and the paper. □ 
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A Appendix: Classical and semiclassical pseudodifferential 
calculi 

A.l Semiclassical calculus 

Our proof of the L 2 estimate in the small frequency region requires the use of classical and 
semiclassical pseudodifferential operators with finite regularity in x'. Here we summarize 
the needed properties of those calculi, referring the reader to the Appendix of [GMWZ2] for 
all the proofs. We are not able to use paradifferential operators (which might have allowed 
us to assume much less regularity in x'), because the process of paralinearization introduces 
0(\U\ip) errors at a stage when they are too big to be absorbed by the left side of our 
degenerate I? estimate. 
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Notation A.l. 1. Let Q 1 = (CoiC") £ denote variables dual to the tangential variables 
x' = (xq,x"), and set ( = where we always take 7 > 1. Set (Q = ylCI 2 = 7p 

and, with slight abuse, ((') = 1| 2 - 

2. For e > /et /3 = (/3',y) G R d x R + fresp. G R d j denote a placeholder for e( 
(resp. e('). 

3. We will ignore powers of 2tt in all formulas involving pseudodifferential operators 
and Fourier transforms. 

4- On H s (R d ) define the norms \u\ Sjl = |(C) s ^|l 2 - 

5. The notation 

for a family of linear operators mapping one function space into another means that the 
operator norm is uniformly bounded with respect to 6,7 for < e < 1 and 7 > 1. For a 
particular s G R we say T £i7 is of order k on H s if 

T ea : H s {R d ) -»• H s ~ k (R d ). (A.l) 

When the domain and target spaces of T are clear from the context, we'll write simply \T\ 
for the operator norm. 

6. We will sometimes denote spaces like C M (R d ,,C°°(R d x R+)) by C M {x',C 00 (l3)) 
when the domains of the variables involved are clear. 

Remark A. 2. Our pseudodifferential operators are defined by symbols with finite regularity 
in x'. Such an operator is generally of order k on H s only for s in a proper sub interval of 



The semiclassical operators are built from "symbols" in the set 

S M = {p(x',p) G C M (Ri,,C oc (R d xR + )) : 

p is independent of x for \x \ large and sup \d^,dp,p(x' , (3)\ < C v }. 

(A.2) 

Let iSqo = (1m<Sm- Define symbol norms 

\p\m,k = sup sup sup \d%,dg,p(x' , /3)\. (A.3) 

\fM\<M\u\<K (x',/3) 

To each p(x f , (3) G Sm we associate the operator defined by 

p(x', eD)u = J e ix 't'p(x', eCMC'R'- (A.4) 

Proposition A.3. If p G Sm and M > d+1 then 

p{x',eD) : L 2 (R d ) -»• L 2 (R d ). 

Definition A.4. A family of linear operators r £i7 is said to be of order e k if r £i7 = e k TZ £jl 
where 

K en : L 2 (R d ) -»• L 2 (R d ). 
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Proposition A. 5 (Products). Suppose p G Sm± and Q S Sm 2 > where M\ > d + 1 and 
M 2 > Mi + (d + 1) + k + 1 for some k>l. Set 

t(x',P)= £ ^ al d%,p(x',P)D«,q(x',P). (A.5) 

|a|<fc-l 

Then t(x' , (3) G Smi and 

A = p(x', eD)q(x', eD) = t(x' , eD) + r e , 7 , (A.6) 
where r £i7 is o/ order e fc . Precisely, r £i7 = e fc T, where 

\T\ < C\p\d+i,k\d x 'q\M 2 -i,o- 

Proposition A.6 (Adjoints). Suppose p£Sm, where M > (d+ 1) + k + 1, for some k > 1. 
Se* 

t ( x '^)= E ^ H ^V(*',/3). (A.7) 

|a|<fc-l 

T/ien i G and 

p(a/, eD)* = t(x' , eD) + r £)7 , 
where r £j7 is o/ order e fc . VFe /taue r e)7 = e fc T ; where 

\T\ < C\d x ,p\ M -i,k. 

A. 2 Classical calculus 

For m G R define the classical symbol classes 

C™ = {p( x \ Q G C M (]Rf,, C°°(M d x {7 > 1}) : p is independent of x' 

for \x'\ large and sup |d^,p(x, C)| < C , I/ <C) m " 1 " 1 }, (A ' 8) 

and set = C\mC^. Define associated symbol norms 

\ P \m,k= sup sup sup|^^p(x',C)KC) M " m - (A.9) 

\/j,\<M\u\<K(x',t) 

To an element p(x',Q 6 C]J} we associate the classical operator 

p(x',D)u = J e ix '('p{x'X)u{C)d£. (A.IO) 

Proposition A.7 (Classical products). Suppose 

p(x',C)GQ andq(x',()eC™l, 

where M x > d + 1 and M 2 > 2(d + 1) + |mi| + 3. Sei i(x',C) = p(^', CM^', 0- I^en 
t G C^ +m2 and 

p(x / ,D)g(a/,D)=t(x / ,D) + r, (A.ll) 
where r is of order m\ + m.2 — 1 . We have 

\r\ < C\p\ d+h i\d x >q\ M2 -i,o- 
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Proposition A. 8 (Classical adjoints). Suppose 

where M > (d + 1) + \m\ +3. Set t(x', C) = p*(x', C). T/ien 

=t(a/,D) + r ei7 , (A.12) 
where r is of order m — 1 and |r| < CI^/pIm— 1,1- 

A. 3 Garding inequality 

Notation A. 9. i. Lei (n, n) denote the L 2 pairing, which can be extended as the duality 
pairing on H s x H~ s . 

2. For a matrix a (symbol or operator) set 3fta = . 

The following Garding inequality is used in the proof of the 1? estimate to obtain bounds 
from below both in the interior and on the boundary. 

Proposition A. 10 (Garding inequality). Consider n x n matrix symbols a G Cj^, w G 
C° M2 , where M x > 2(d + 1) + max(f ,m) + 2 + [|f |] and M 2 > 2{d + 1) + m + 2 + [|f |]. 
Suppose there is a scalar symbol \ £ and c > snc/i i/iai x 2y; = w an d 

&a(x', C) > c(C) m on supp x- (A.13) 

Let A = a(x' , D) and W = w(x', D). Then there exists C > such that for all u G H~ 

< 5ft(^Wn, VFn) + C|n||_ li7 . (A.14) 

The constant C depends on symbol norms of a, w, and \. 

B Appendix: An existence result for systems with pseudod- 
ifferential boundary conditions 

In this appendix we discuss how to solve variable-coefficient hyperbolic initial boundary- 
value problems with pseudodifferential boundary conditions, such as arise from applying 
the procedures described in Section [1.51 to the quasilinear case. The problem (|1.6U|) is also 
of this type, but with constant coefficients. 

Consider a noncharacteristic hyperbolic boundary-value problem on the half-space {x = 
(x',x d ) = (t,x",x d ) : x d > 0}: 

u Xd -G(x,D')u = f€e' t L 2 (x), 

e^T(x\D')e-^u{x'^) = g G e' l l?{x'), { 
where G is a differential operator with symbol 

d-l 

G(x,iT,i v ) := -A^\x)(iT + J2iVjM x ))> ( B - 2 ) 
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derived from the hyperbolic equation u t + Aj(x)d Xj u = Ajf. Here the Aj(x) G C°° 

are N x N matrices and constant outside some compact set. We suppose A<i{x) is invertible. 
Assume also that the eigenvalues A ? -(x, £) of Sj=i are C°° , real, semisimple, and of 

constant multiplicity for all (x, The boundary operator is a classical pseudodifferential 
operator of degree zero associated to a C°°, p x N matrix symbol T(x',t, 7, 77) using the 
quantization of section IA.2I 

Because of the pseudodifferential boundary conditions this problem is not covered by 
the standard theory presented, for example, in |CPl IBS] . To state our result, we first define 
D(x', t,7, rj) to be the Lopatinski determinant 

r, 7, 77) =det(kerr(x', r, 7,??), E + (-G(x', 0,r, 7 ,7/))) , (B.3) 

where the spaces appearing in the determinant are defined just by freezing x' . Here, follow- 
ing [Met4[ pp. 135-136], we mean the determinant obtained by substituting orthonormal 
bases of the spaces involved, so that the modulus of the determinant is independent of the 
choice of basis. 

Proposition B.l. With the above assumptions on (IB.ip . suppose also that the uniform 
Lopatinski condition holds: 

\D(x , r, 7, 77)1 > C > 0, for C independent of(x,T, rj) and 7 > 0. (B.4) 

Moreover, suppose there hold also the uniform bounds 

|r|,|rt|<c, (b.5) 

where := T* (IT*) -1 is the pseudo-inverse of V. Then, there exists 70 such that for 
7 > 70 the problem (jB.ip has a unique solution satisfying 

7l< 7 + « 7 < C + (g)l 7 ] , (B.6) 

where the constant C > is independent of 7. 

Proof. 1. The proof of the a priori estimate (|B.6P may be carried out using Kreiss sym- 
metrizers as in \CP\ IBS] , since 9 Xd — G(x,D') may be conjugated to block structure and 
the uniform Lopatinski condition holds. See also [Met 41 pp. 135-136], particularly Lemma 
6.2.4, for details of the Kreiss symmetrizer argument and the role of bounds (|B.5p , 
2. Next we define the approximate adjoint problem 

u Xd + G(x, D'fu = f e e-^L 2 (x), 
e-^t{x',D')e^ t u{x',Q) = ~g^e-^L 2 {x') 

where G(x, D')* is the formal adjoint of G(x, D') and the (N — p) x N matrix symbol T 
is chosen so that Kerf = (Kerr) -1 - for each (x f , r, 7, rj). Without loss of generality, we may 



19 Matrices Aj(x) of finite regularity can be treated by similar arguments if one uses paradifferential 
operators. Our assumption on the Xj, which implies that the operator dd — G(x,D') in (|B.1|) can be 
conjugated to block structure in the sense of |MZ2| . can be weakened as discussed in MZ2 . 
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take both T and T to have orthonormal columns, so that the symbols satisfy Id = T*r+r*r. 
Let us now write T(x',D') = T D ,, G(x,D')* = (G D >)*, etc.. Since 

(r D ,yr D , = (r*r) D , - r d , and 

(f D ,yf D , = (f*f) D ,-R D ,, ' 8) 

where Rd 1 and Rpi are operators of order —1, we have 

Id = (T D ,)*T DI + Rjy + (f D ,)*f D , + Rd,. (B.9) 

3. Observe that for 7 large we can invert the N x N matrix operator ( ^ D J exactly 



on I? by first using the calculus to construct an approximate inverse, and then using a 
Neumann series. Thus we can solve the operator equation 

+ Rjy = ((S D 'T (Sd'T) f = (Sd>)*Td> + (S D >yf D > (B.10) 

for (Sd>)* and (Sd 1 )*- This gives 

id = (r D > + s D >yr D > + (f D < + s D ,yr D ,. (b.h) 

The operators Sd' and Sd 1 have norms of size 0{^~ 1 ) on I? . 
4. Let us define 

T D , n ■= e 7 *r D , e - 7 * and T D / ) _ 7 := e^Vjye^. (B.12) 
Integrating by parts we obtain the duality relation 

((d Xd -G D ,)u,v) - (u,(-d Xd - (G D ,y)v) = (u,v) = (e 7t n,e- 7 ^) 

((r£>/ )7 + Sd'^u^d',-^) + ((ro' i7 + s , D/ i7 )u,r £) / _ 7 u) 
G c c °°(I 

on (d Xd + {G D >)*)E by 



(B.13) 

Setting £ = C C °°(I + ) : r D / _ 7 u = 0} and recalling ([BTf]) . we define a linear functional 



i ((-d Xd - {Gjy)*)v) := (/, v) - (g, Tjy^v). (B.14) 

One checks that the forward Lopatinski condition (|B.3|) implies that the adjoint problem 
(IB.7P satisfies the backward uniform Lopatinski condition (see [CP| ). and thus 



7M0-7 + M°))o -7 < 1 for v e £?, (B.15) 



where |it|o,- 7 := \e yt u\ L 2. Using 

\{g,T D ,^v)\ <C(g)o„(v) ,- y , (B.16) 
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a similar estimate for (f,v), and (|B.15p . one obtains readily that i is continuous on (d x , + 
(Gj)r)*)E for the topology induced by e _7 *L 2 . The Riesz theorem then yields u G e 7 *L 2 
such that 

(u,(-d Xd - (Gjy)» = (f,v) - {g,T D ,^v) for all v G E. (B.17) 

Taking u G C%°{xd > 0} we obtain from (|B.17|) that d Xd u — Gdu = f. In view of 
pT[%pr[7j) we obtain 

((ri)',7 + Sd',~()u — 9, = for all v £ E. (B.18) 

Since the matrix symbol T is surjective when restricted to kerT, it follows (by density) that 

(Td',7 + Sd'^)u = g. (B.19) 

5. Given (f,g) G Y 1 := e 7 *L 2 (x) x e 7 *L 2 (x'), we have found u G e 7 *L 2 satisfying the 
nearby problem 



u 



X d 



G(x,D')u = f G e^L 2 (x), 



{Y DI)1 + 5 D , i7 )n(x',0) = g G e 7i L 2 (x'). 
Applying the estimate (|B.6P and treating Sd>^u{x' , 0) as an absorbable error, we have 

l/|o, 7 



(B.20) 



(n)o, 7 < C ( ™ + (5)0,7 ) , (B.21) 



C f|/|0, 7 



and thus 

<fik', 7 «>o, 7 < ^ + (<7}o, 7 J • (B.22) 
6. Let T~ l denote the inverse we have just constructed for the operator 

T a u := L ^? " Gd '^ ) , and set Tu := ( ®- " Gd '>) . (B.23) 

\( L D',-y + JD',y)U\x d =0/ \ >- D' ,-yU\x d =0 J 



We have T = T a + _ , so TIT 1 = I + _ T~ L := I + M T The estimate 

y — JD^-yJ \—JD',t/ 

(|B.22p implies that M 7 has norm < 1 on 7 7 for 7 large, so we can invert I + M 7 on Ky by 
a Neumann series, and thereby invert T. 

□ 



C The example of Rao 

We now consider the example left untreated in the analysis of (jl.ip in [R], lying in the 
case (ii) described in Section 11.41 Consider perturbations about a constant boundary-layer 
solution of (11. lj) in dimension d = 2, with nonvanishing tangential velocity 

(C.l) 
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and normal velocity v noncharacteristic, i.e., ^ v,v ± c, where c := yPp + ^p~f" is sound 

speed. Here, we are assuming a polytropic temperature law e = c v T (used but not stated 
in Chapter 5 of [R]) and an unspecified pressure law p = p(p,T), with p p + > 0. 

Working with variables U := (p, u, v, T) T , we find, following [PQ, that the equations (jl.ip 
may be expressed in quasilinear form as 



A d t U + ^ Ajd Xj U = eJ2 d Xj (B jk d k U), 

3,k 



(C.2) 



with 
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Pt/p 



(C.4) 



and 



B 



■jk 



(C.5) 



(0,0,0)' 

o p jh 

where f3jk is elliptic in the sense that the eigenvalues of Y?,jk PjkCjCk have real part > c|^| 2 , 
c> 0, for all ^ e R d . 

Note that the parabolic terms involving Bjk are of the more general form treated in 
[GMWZ5, GMWZ61IR] and not the Laplacian form to which we have restricted for simplicity 
in (ll.3ip . In particular, they are degenerate parabolic, no diffusion being present in the 
equation for density p, whose principal part pt + (u, v) T ■ V ' x p is a hyperbolic transport 
equation. Though it does not change the theory in any essential way, this will require a bit 
of discussion as we go along. 

To start with, and most significantly, the fact that the p equation is hyperbolic means 
(see (SHU EU |U IGMWZ51 IGMWZ6] ) that a boundary condition is needed for p only when 
the convection v in the normal direction is inward on the boundary, i.e., v > 0. Boundary 
conditions on the "parabolic variables" (u, v, T) must always be imposed, and may be 
Neumann or Dirichlet type. Here, we will impose Dirichlet conditions on p (when needed), 
u, v, and Neumann, or "insulative" conditions ddT = on T, following (one version of) 
engineering practice as described in the introduction. 

Assume, now, that we are in the incoming supersonic case in the notation of [GMWZ5, 
|R] . the case left open in the treatment of [Rj, that 



< c < v. 



(C.6) 



By the previous discussion, we require a full set of boundary conditions for the hyperbolic- 
parabolic problem, including a condition on p. We thus take boundary conditions 



(p,u,v) = (51,52,53) and d X2 T = = 54 at x d = 0. 



(C.7) 
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The first-order hyperbolic part of (|C.2p . comprising the Euler equations, may be written as 



U t + J2Ajd Xj U = Q, (C.8) 

3 

where Aj := Aq 1 Aj. Recalling the standard computation that the characteristics of A2 are 
v — c, v, and v + c (with multiplicity two), we find, consulting (lC.6p . that this hyperbolic 
problem is totally incoming, i.e., all eigenvalues of A2 are strictly positive. 

Though we shall not carry it out here, it is not difficult to see using the methods of 
[GMWZ5, GMWZ6J that in this situation there can exist no small-amplitude boundary 
layers other than the trivial, constant layer, for the simple reason that any rest state sat- 
isfying (IC.6P must, by the dimensional counting arguments of [GMWZ5I, be a repellor for 
the standing-wave ODE, so cannot be the limit as x<i — > +00 of a nonconstant standing 
wave (boundary layer). Thus, the same derivation as in Section 11.41 case of a for- 
mal boundary-layer expansion applies, yielding an outer problem with D = 3 Dirichlet 
conditions and N = 1 Neumann boundary conditions, the same ones imposed on the full 
hyperbolic-parabolic problem. 

At this point, having derived an outer problem, we can forget its hyperbolic-parabolic 
origins and analyze it as in Section 11.41 As noted in Lemma II .11( d). we have that the 
outer problem is weakly Lopatinski stable, and solvable for the constant-coefficient problem, 

provided 7-2) or, equivalently, f^ 1 ^? is full rank, where 

\ L 2 A 2 J \ ±2 J 

/l 0\ 

Ti := 1 , T 2 := (0 l) . (C.9) 
\0 1 0/ 

are the boundary matrices corresponding to the boundary conditions described above. 
Computing, we have 

fT 1 A^ 1 A 2 \ _ ( M * N 



1 \ / v p\ 
where M = [ — u/p 1/p II v I is invertible provided that 
-v/p 1/p) \ Pp /p v) 

(v p\ 
v J = v{v 2 -p p ) = v 2 - c 2 +ppr/p 2 , 
Pp/p v) 

or 

v 2 — c 2 + ppx I p 2 ■ (Weak Lop) 

By the supersonic condition (lC.6p . this is evidently true under the standard assumption 

p T > 0, (CIO) 

satisfied in most typical applications, in particular, for an ideal gas pressure law p = RpT, 
where R > is the universal gas constant. However, in general (jC.lOp . hence apparently 
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also (Weak Lop), need not be true, for example for nonstandard materials/phases such as 
ice, for which pressure decreases with temperature. 

Continuing, let us check not only reality but also semisimplicity /constant multiplicity 
of the characteristic roots of the boundary problem 

(Ida 0)\ /(Ida 0)\ 

i / )dtv + vn[^ ._i / )div = 



derived above, or, equivalently, eigenvalues of 

" 171 {r 2 A^A ) [tiA^aJ -- i7] {* (T^A.Tlr^A^A.Tl) ' ^ U > 

where, by our previous calculations, T 2 A2~ l AqTT} ^ 0. By the lower block triangular struc- 
ture of (jC.lip . this is implied by the nonvanishing property T 2 A 2 ^iT^ / 0. Noting that 

/ \ 



1 2 



Pp/P 



V u ) 



while, by Kramer's rule, 



V 2 A 2 X = (detA 2 ) 1 (vp T p/p 2 c v -v 2 p/pc v (v 2 - p p )v) , 



we find that I^Ag 1 Air^ = uv(v 2 —p p )/ det A 2 is nonvanishing if ( |Weak Lop ) holds, under 



the nonvanishing tangential velocity assumption (jC.ip . The condition of constant multi- 
plicity is trivially satisfied, since rj is dimension one. 

On the other hand, when (IC.lj) is violated, we have, computing, 

T 2 A 2 ~ 1 A 1 = (detA 2 y 1 (* v 2 (p(v 2 - Pp )-p/c v ) 0) , 

which in general does not vanish, and so the matrix on the righthand side of (jC.lip contains 
a Jordan block, violating semisimplicity. (Specifically, it has all zero eigenvalues, by lower 
triangular form and vanishing on the diagonal, but is not identically zero.) This shows in 
passing that semisimplicity (hence maximal bounds) can fail for method two when the weak 
Lopatinski condition is satisfied. 



Thus, when (Weak Lop) holds, in particular under the standard assumption (IC.lOj) . we 
have also weak hyperbolicity together with semisimplicity/constant multiplicity of charac- 
teristics of the boundary problem, so obtain by the theory of Section Tl .61 (together with the 
observation above that boundary layers are absent at all orders in this case) existence of 
approximate solutions to all orders. However, the question of convergence is still open up 
to now, for lack of associated hyperbolic-parabolic estimates. 

Remark C.l. Example 11.21 is closely related, and gives a similar conclusion, as does any 
problem with a single Neumann condition. 

D Extension of the second approach to the non-totally in- 
coming case 

We return now to the reduced hyperbolic problem in its original form (jl.55H in the 
general case of mixed-type boundary conditions, assuming that A^ is nonsingular and that 



56 



the Aj are symmetric. We also assume that the operator L can be conjugated to block 
structure in the sense of [MZ2J. 

We split the problem (jl,55p into two parts: a problem with homogeneous Neumann 
boundary conditions 

Lv = f, 7T + (A d )d d v\x d =o = 0, (D.l) 
and a problem with homogeneous forcing and mixed boundary conditions 

Lw = 0, Tiw\ Xd=0 = gi -Tiv\ Xd=0 := hi, t 2 d d w\ Xd=0 = h - f zdd,v\ Xd =o := h 2 (D.2) 

where uq := v + w. In t < we have v = and w = 0. 

The first problem may be solved as in Section Tl. II To solve the second problem we take 
the Laplace-Fourier transform, substitute the usual boundary symbol for d d , and multiply 
r"l by (7 + ir) to obtain 

(7 + \ _ /( 7 + i T )h 



12 

Note that to is a decaying solution of Lw = if and only if w(0) € K + (A' d 1 (j + it + 
2Zj=l iVjAj) j^l- Thus, the problem Lit; = on x d > with boundary conditions (ID.3j) is 
equivalent to the problem on the boundary with enlarged boundary conditions 

r (7,r,7 ? )(7 + ir) \ /0\ 

ri(7 + *r) u>(0) = Ui , (D.4) 

-f 2 a- i ( 7 + + e •=! ^ao7 w 

where 1^0(7, r, 77) is a matrix whose rows are orthogonal to K+^A^ 1 (7 + it + X^j=i ^Vj-^j, 
But, this is equivalent to the Cauchy problem on the boundary 

To(7,t,??)\ d-i I \ /0\ 

Ti + y~M J to^ = \h x J where /i 3 = -h 2 . (D.5) 

IW / i=i KhAfAj \h 3 J 

Definition D.l. Let = {(^7,77) : |r, 7,7y| = 1,7 > 0}. Parallel to Definition \1.14\ we 
we say that (|D.5|) is evolutionary at (ro,7o,^o) G 5"+ */ 

^r (7o,r ,7?o) N 

,4,, - i ri 

is invertible. We say (|D.5j) is weakly hyperbolic at (70,70,770) G 5"^ «/, in addition, the 
first- order system 

d-i /ro(7o,To,r/ )\ _1 / \ /r (7o,7-o,r?o)\ _1 /0\ 

(7 + »t)u> + £;*7& ri «> = r i M (D - 6) 

i=i V f^- 1 / Vfwv \ fw 7 VW 



20 Here and in solving ()D.1|I we use the assumption of constant coefficients. 

21 By the block structure assumption this space extends continuously to 7 = [MZ3| 
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has pure imaginary characteristics, defined as eigenvalues of the homogeneous degree one 
symbol A(r};~fo, to, rjo) := Ylj=i ^j-^jdo, T o, Vo) with homogeneous degree zero coefficients 

/r (7o,T ,77o)\ _1 / \ / 

A(7o,ro,ryo) := ^•q 1 (7o,t ,t?o)A' = r i > A := 

V ^A d l J \Y 2 A~ l Aj \t2A^A jj 

Unlike the situation of the totally incoming case, the evolutionarity and weak hyperbol- 
icity conditions must be checked separately, and do not follow simply by the weak Lopatinski 
condition. However, we have: 

Proposition D.2. System (1D.5|) is evolutionary and weakly hyperbolic at all (70,70,770) £ 
only if the original problem (|1.55|) satisfies the weak Lopatinski condition. 

Proof. If the weak Lopatinski condition fails for (|1.55|) . then for some 7 > and r, 77, there 
exists w such that ro(7, t, rj)w = 0, Y\w = 0, and f^-Aj 1 ^ + ir + X^/=i ^Vj Aj)w = 0, 
whence (|D.4j) and (1D.5|) hold with /ii = 0, /13 = 0, as, by the evolutionarity assumption, 
does (|D.6p . But, this, by inspection, means that 7 + ir is an eigenvalue of .4.(77; 7, r, 77) that 
is not purely imaginary, contradicting weak hyperbolicity at (7, r, rj). □ 

When the evolutionarity and weak hyperbolicity conditions do hold, we have a situation 
analogous to that of Section 11.61 but for the generalized first-order Cauchy problem 

d-l 

d t w + J2 A 3 d x j w = G, (D.7) 

_ . /r (7,r,7 ? )\- 1 /o 

on the boundary Xd = 0, where AjW := Aj(j,T,r))w and G := Ti I hi 

\ 1-2.1,/ ; / W 

As in the standard case, evolutionarity plus weak hyperbolicity is not sufficient to give 
well-posedness of the problem ID .71 but requires some additional structural assumptions. 
The following sufficient conditions are often applicable. 

Lemma D.3. Assuming evolutionarity at all (10,70,770) £ S + , a sufficient condition for 
well-posedness of (|D.7ji 22 l is that the eigenvalues ofA(n; 70, to, 770) be semisimple, pure imag- 
inary, and of constant multiplicity for all (70,70,770) 6 S + , 77 ^ 0. In in this case we have 
the uniform resolvent estimate 

d-l 

|(7 + ir + ^^ i (7,r,77))- 1 | < C/ 7 , 7 > 0. (D.8) 
3=1 



22 By well-posedness of (|D.7[) we mean that there exists 70 such that for 7 > 70, if G G e 7 *L 2 , then there 
is a solution w £ e Jt L 2 . 
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Proof. By taking the Laplace- Fourier transform of the problem (|D.7|) . we see that well- 
posedness follows from the estimate (|D.8p , We prove (|D.8|) by first proving the family of 
estimates 

d-l 

|(7 + ir + ^^ j ( 7 o,To,%)) _1 | < C/ 7 , 7 >0 (D.9) 

3=1 

for the problems frozen at (70, To, 770) E S + . For a given (70,70,770) this estimate can be 
obtained by conjugating the frozen system to Jordan form and using semisimplicity. Using 
compactness of the unit half-sphere and continuity of the Aj , we conclude that the estimate 

is uniform for all (ro,7o,7?o) £ S+- Restricting to the diagonal (70, To, 770) = (7, t, f?)/|7, t, 77 1 
and using degree zero homogeneity of Aj , we obtain (1D.8|1 . □ 

The conditions of Lemma [D. 31 which refer to a large family of frozen systems, seem hard 
to check and possibly over-restrictive; however, as examples below show, the conditions are 
sometimes satisfied. Since the matrix A has a guaranteed large block of semisimple zero 
eigenvalues associated with the boundary conditions To and T± , in simple cases verification 
amounts to checking nonvanishing of a single remaining eigenvalue (see the proof of Propo- 
sition [1X5]) . For systems of size N = 3, for example, with 2 Neumann conditions, 1 Dirichlet 
condition, and 2 incoming modes, the matrix A has a two-dimensional kernel, so we are in 
the situation described. Likewise, for systems with < 2 Neumann boundary conditions, if 
there is even a single outgoing mode, then we are in this case. Of course, one might always 
attempt to establish (|D.8|) by direct matrix inversion, as well. 

When we do have conditions sufficient to obtain well-posedness of (|D.5p . hence bounds 
(|D,8j) on 1^1^=0, we may then from the trace information on w \ Xd =o obtain interior bounds 
on w using any Lopatinski stable Dirichlet boundary condition, for example a maximally 
dissipative one. Putting the estimates for v and w together, we obtain a Kreiss-type estimate 
with losses for the solution uq, similarly as in the totally incoming case. 

Remark D.4. Even when the semisimplicity and constant multiplicity assumptions of 
Lemma ID. 31 are not satisfied, it may still be possible to implement method two in some 
problems. This would require refined estimates for the Cauchy problem that to our knowl- 
edge have not yet been carried out. 

D.0.1 A sharp condition in a special case 

Let the number of Neumann conditions be one greater than the number of incoming hyper- 
bolic modes, i.e., the number of reduced Neumann conditions be one. 

Proposition D.5. Let the number of reduced Neumann conditions be one. Then, assuming 
evolutionarity at all (r, 7, 77) E S + , a necessary and sufficient condition for well-posedness of 
(|D.7j) (i.e., for the estimate (ID. 8ft ) is that the eigenvalues of A{n; 0, 0, rf) be semisimple, pure 
imaginary, and of constant multiplicity with respect to r\, for all \r]\ = 1, or, equivalently, 
the scalar condition Ylj^d^Vj^j f or V 7^ 0, where ctj := (AjA$ 1 )atat(0, 0, rf). 

Proof. Conjugating by the uniformly invertible Ao, we have that (|D.8P is equivalent to 

K(7 + ^)+E j ^^Mi(7,T,7 ? )^- 1 )- 1 |<C/7for7>0, (D.10) 
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or, computing explicitly in block-diagonal form, 



C/l> 



7 + %T 



>.],,,/ *VjPj 7 + + Y j:d irjj 



a 



-i 



7+JT 







(D.ll) 



\ 7+ir (7 + »^ + Ejyi«w) / 

where the upper blocks consist of N — 1 rows, the lower blocks of 1 row, and 

Pi = (/3},...,/3f- 1 ), $( 7 ,t,7/) := (^VK<(7,t,»7) 



(D.12) 



are real and scalar. This holds for Y),t\ « \rj\ if and only if Ej^d ^ r li ot j{^-> 0> 7 ?) 7^ f° r 
r/ 0, or, equivalently, | Ej^d 7 ?.? ?^' ^/Ol > M/C > 0- For I 7 ?! — C|^, x| , it holds always, 
in the 1 — 1, 1 — 2, and 2 — 2 blocks by inspection using the fact that 7 + it + Ylj^d^i 01 ! 

7+ir 



is scalar, and in the 2 — 1 block by 



<CH/| 7 ,r| <C 2 . 



□ 



Remark D.6. Though the proof of Proposition ID. 51 relied heavily on the scalar nature 
of block 7 + it + Ylij^difljCXji the same argument shows in the general case that, assuming 
evolutionarity, necessary and sufficient conditions for (ID.Sjl are well-posedness of the reduced 
symbol in the 2 — 2 block, 



|(7 + *r + . .,,"/,'!;; x | < C/7 for 7 > 0, 



(D.13) 



and noncharacteristicity at (7, r) = (0,0), det(^ y d ir]jOtj(0, 0, 77) ) 7^ 0, for 7/ 7^ 0. Though 
still nontrivial, the formulation (|D.13|) represents a substantial reduction in order. 



D.l Some illustrative examples 

We now give some examples indicating the range of possible behaviors. Before presenting 
these, we first prove two theoretical results that serve to frame the discussion. 

Lemma D.7. Assuming that Aq is invertible for (r, 7, 77) E S+ and that 



M 



F 2 A 2 (7 + *r + Ej#d 



is full rank for 7 > and I7, r, 77 1 = ljff|, (D.14) 



i/ie uniform Lopatinski condition for the rescaled boundary condition of method one (|1.62|) 
may be expressed, equivalently, in the convenient form 

^ + iT y a nkt 2 -N det ^ + iT + j2 j ^ d i Vj A j (j,T,rj))\ >5 >0 for all 7 > 0, |7,r,7?| = 1. 

(D.15) 



23 For pure Neumann boundary conditions, and non-totally incoming hyperbolic characteristics, this holds 
generically for dimension d = 2 in the sense that it is true for all choices of boundary matrices V2 except for 
a measure zero set, and fails generically for d > 3 in the sense that it is false away from a measure zero set 
of matrix entries Ai, . . . , Ad (precisely, those for which Ai, . . . , Ad-\ all share a common eigenvector). 
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Proof. Under (|D.14|) . \M\ and its pseudo-inverse := M*(MM*)~ 1 are uniformly bounded 
on 7 > 0, I7, r, 77 1 = 1, by continuity, as is Tq (which has orthonormal rows) and its pseudo- 
inverse. Thus, (see [Met4, pp. 135-136]) the uniform Lopatinski condition, or transversality 

/ r N 



of the kernels of M and To, is equivalent to | det 



M 



> 5q > 0, whence equivalence of 



(|D.15|) follows by uniform boundedness of det A 1 (a consequence of continuity and evolu- 
tionarity) and 



del 1 Ti^ + ir) 

k^A- 1 ^ + it + ^2 j ^ d ir jj A j 



(7 + it) 



rankAf-r 2 



det 



r 

,1/ 



□ 



Corollary D.l. Under the assumptions of Lemma \D. 7\ the uniform Lopatinski condition 
is sufficient (but not necessary; see the examples below) for resolvent estimate (|D.8j) . 

Proof. Conjugating by the uniformly invertible Aq, we have that (|D,8|) is equivalent to 



\((l + ir) + Ej^d iViMli r > V)A Y 1 ] < Ch for 7 > 0, 
or, computing explicitly in block-diagonal form, 

-1 / -L- ^ 



Ch> 



7 ~\~ VT 

l,T,r]\/3 \j,T,i]\a + 



7+ir 
7+ir 



(D.16) 



(D.17) 



where a+ := 1+lT+ ^p t i ^ %ajr,i and f3 := 1 a j an d Pj homogeneous degree zero in 

(7, t, 77), defined as in (|D.12p . By (|D.15|) . Uniform Lopatinski is equivalent to |detaq_| > 
5q > 0, whence, by boundedness of a + (j,r, 77), ja^ 1 ) is uniformly bounded. This, along 
with boundedness of /3, verifies (|D.17p . □ 

Example D.2. Consider the system ut + A\u Xl + A2U X2 — eA x u = f, u S M 3 , with two in- 
coming hyperbolic modes, two Neumann conditions T2u\ X2 = g<i, and one Dirichlet condition 
ri7i| X2=0 = 92, given by 





1 





















1. *H 





1 




Vo 





I) 




Vo 








and r 



1 



We have evidently T2 



(D.18) 

(0 a P). More, by 



1 * * N 
v a 

the decoupled form of A\, A2, and the fact that H + (A2)u = is a maximally dissipative 
Dirichlet condition, we find without calculation that K+(A2 + ir + 777^2)) is spanned by 
(0,0, 1) T and (*,1,0) T , so that 



Aq :-- 



To 
S2A2 



'1 5(7, r, 77) 0^ 
10 
,0 a (3 1 



A! :-- 









(D.19) 



61 



Computing, we find that evolutionarity, Aq invertible is satisfied when f3 7^ 0, with 



(i -6(rf,T,Tj) o \ / o o o\ 

Aq 1 = 1 , A l Ai = 0, (D.20) 

\0 -a/13 1//3J \a//3 1/ 

in which case .Aq 1 ^! is independent 0/(7, r, rj), with real, semisimple constant-multiplicity 
characteristics 0, 0, 1. Thus, by Lemma \D.3\ method two works in this case. 

On the other hand, by (|D.15P and the easily-checked (|D.14p for a / 0, the uniform 
Lopatinski condition for the rescaled boundary condition of method one is equivalent to 

h + ir \ 

< £0 < |(7 + ir)~ 2 det 7 + ir | = I7 + %{t + rj)\ 

\irja/(3 j + ir + irj/ 

for 7 > 0, |7,r, 77 1 = 1. Tftis clearly fails for 7 = 0, 77 = — r. Thus, method two can apply 
even when uniform Lopatinski fails. 

Example D.3. Consider again Example \l.ll of the first-order wave equation with drift a, 

Al ={l oj' ^ 2= ( -1 + aJ' H<1 ' 
with full Neumann boundary conditions, so that T% = (l 0), To = (£(7, t, 77), 1). Then, 
f 2^-2 1 = dfe 0) and = (0 ^) , so tftat A = (j_ JY = (° _?_Y 

4-1 (1 + a) \ . Jt ._, . A) 1 \ . . , , x . ± , 

A = I ^ J(i + ck) / ' a S = lo 5 J lS a ~ m "' e P en dent, with real, semi- 

simple eigenvalues of constant multiplicity whenever 5 7^ ; or, equivalently, 77 7^ 0. Applying 
Lemma \D.3\ we thus find that the associated Cauchy problem of method two is well-posed 
with standard hyperbolic estimate (|D.8p . Thus, again, method two succeeds despite failure 
of the uniform Lopatinski condition observed in Example \l.l\ for method one. 

Example D.4. Substituting in Example \D.3\ the value A± = ^ ^\ , 6 7^ 0, we find that 

AiAq 1 = ( 1 „ ~\, and so semisimplicity fails, in general, for r/ 7^ and \ j, r\ 7^ 0. 
\T+a v — oj 

Thus, in this case the frozen- coefficient algorithm proposed for method two fails, even though 
by Proposition ID, 51 the problem is well-posed. 



D.2 Questions and comparison to first approach 

We have described two methods for solving the reduced hyperbolic outer problem ob- 
tained by the derivation of Section 11.41 which appear to give slightly different bounds and 
apply in slightly different scenarios. It would be interesting to further clarify the relation 
between the two methods. It is not clear that our technique of obtaining estimates through 
an enlarged family of frozen-in frequency coefficients -4,o(70i tq,t]q) will always produce opti- 
mum results for the problem of method two in problems of interest. On the other hand, for 
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situations with mixed incoming/outgoing modes, the first method requires that the uniform 
Lopatinski condition be satisfied in order to obtain good bounds, a scenario that might not 
occur even when the conditions for method two do apply. 
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